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We perform first-principles calculations of the structural, electronic, vibrational and magnetic 
properties of the C4SN12 azafullerene and Ceo- Full geometrical optimization shows that C4SN12 is 
characterized by several distinguishing features: only one nitrogen atom per pentagon, two nitrogen 
atoms preferentially sitting in one hexagon, Se symmetry, 6 unique nitrogen-carbon and 9 unique 
carbon-carbon bond lengths. The Mulliken charge analysis indicates that the doped nitrogen atoms 
in C48N12 exist as electron acceptors and three- fourths of the carbon atoms as electron donors. Elec- 
tronic structure calculations of C48N12 show that the highest occupied molecular orbital (HOMO) 
is a doubly degenerate level of a g symmetry and the lowest unoccupied molecular orbital (LUMO) is 
a nondegenerate level of a u symmetry. The calculated binding energy per atom and HOMO-LUMO 
energy gap of C48N12 are about 1 eV smaller than those of Ceo- For both C48N12 and C60, the total 
energies calculated with STO-3G, 3-21G and 6-31G basis sets differ from the 6-31G* basis set results 
by about 1.5%, 0.6% and 0.05%, respectively. Because of electron correlations, the HOMO-LUMO 
gap decreases about 5 eV and the binding energy per atom increases about 2 eV. Our vibrational 
frequency analysis predicts that C48N12 has in total 116 vibrational modes: 58 modes are infrared- 
active (29 doubly-degenerate and 29 non-degenerate modes) and 58 modes are Raman-active (29 
doubly-degenerate unpolarized and 29 non-degenerate polarized). It is found that C48N12 exhibits 
eight 13 C and two 15 N NMR (nuclear magnetic resonance) spectral signals. In comparison to iso- 
lated carbon or nitrogen atoms, an enhancement in the dipole polarizability is found due to the 
delocalized -k electrons in C48N12 and Ceo- The average second-order hyperpolarizability of C48N12 
is about 55% larger than that of Ceo- In addition, the effects of basis sets are discussed in detail, 
and the different methods for calculating nuclear magnetic shielding tensors are compared. Our 
detailed study of Ceo reveals the importance of electron correlations and the choice of basis sets in 
the first-principles calculations. Our best-calculated results for Ceo with the B3LYP hybrid density 
functional theory are in excellent agreement with experiment and demonstrate the desired efficiency 
and accuracy needed for obtaining quantitative information on the structural, electronic and vibra- 
tional properties of these molecules. Our results suggest that C48N12 has potential applications as 
semiconductor components, nonlinear optical materials, and possible building blocks for molecular 
electronics and photonic devices. 



I. INTRODUCTION 

Graphite is a stable and abundant solid form of pure 
carbon [1] . In this form, three valence electrons of each 
carbon atom form three strong sp 2 trigonal bonds to 
three nearest neighbors with an equal distance of 0.142 
nm, while the fourth valence electron from different car- 
bon atoms interacts by weak n bonds perpendicular to 
successive sheets with an inter-plane distance of 0.34 nm 
[2] . Diamond is another slightly less stable and less abun- 
dant crystallographic form of pure carbon [1], where each 
carbon atom is covalently bonded to four neighbors via 



sp 3 hybridization at the apexes of a regular tetrahedron. 
In 1985, a fascinating molecule, named Ceo (a truncated 
icosahedron with 20 hexagonal and 12 pentagonal faces, 
and 60 vertices, each of which is at the intersection of two 
hexagonal and one pentagonal faces) was discovered by 
Kroto et ai. [3], and a new form of pure carbon, named 
fullerenes [4], was born. 

Fullerenes can crystallize in a variety of three- 
dimensional structures [5,6], being made from an even 
number of three-coordinated sp 2 carbon atoms that ar- 
range themselves into 12 pentagonal faces and any num- 
ber (> 1) of hexagonal faces [4]. The macroscopic syn- 
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thesis of soot [5] , which contains Ceo and other fullerenes 
in large compounds, plus the straightforward purification 
techniques of the soot which make the pure fullerene ma- 
terials available, have led to extensive studies of fullerenes 
[7-10,12]. 

Doped fullerenes have also attracted a great deal of 
interest due to their remarkable structural, electronic, 
optical and magnetic properties [7-10,12]. For example, 
the doped fullerenes can exhibit large third-order opti- 
cal nonlinearities [9,10] and be ideal candidates as pho- 
tonic devices including all-optical switching, data pro- 
cessing, and eye and sensor protection [9,10]. Another 
example is alkali-doped Cqq crystals, which can become 
superconducting [13,14], for example, at a critical tem- 
perature T c — 30 K [13]. In addition to the endohcdral 
doping (inside fullerenes) and exohedral doping (outside 
fullerenes), there is another type of doping, named sub- 
stitutional doping, where one or more carbon atoms of 
fullerene are substituted by other atoms [7-12], due to the 
unique structural and electronic properties of fullerenes. 
Because boron and nitrogen bracket carbon in the pe- 
riodic table, much attention has been paid to alternate 
boron- and/or nitrogen-doped compounds [7-12]. Over 
the past 10 years, boron and nitrogen atoms have been 
successfully used to replace carbon atoms of Ceo and syn- 
thesize many kinds of hcterofullerenes, C6o- m -nN m B n 
[7-10,12,15-18]. In 1995, a very efficient method of syn- 
thesizing C59N was reported [16]. This method has led to 
a number of detailed studies of the physical and chemical 
properties of C59N [8-10,12,17]. Very recently, Ceo with 
more than one nitrogen atom replacing carbon atoms in 
the cage has been synthesized by Hultman et al. [18], 
and the existence of a novel C48N12 aza-fullerene [18-20] 
was reported. Hence, it would be interesting and use- 
ful to investigate and predict the structural, electronic, 
vibrational and magnetic properties of this aza-fullerene 
by performing detailed first-principles calculations. This 
forms main purpose of the present paper. 

Fullerenes have been challenging molecules for first- 
principles calculations because of their size [21,22]. Re- 
cent advances in ab initio methods and parallel comput- 
ing have brought a substantial improvement in capabili- 
ties for predicting the properties of large molecules. The 
coupled cluster method has been used to predict phenom- 
ena in C20 [24]. Other first-principles methods, which 
are less demanding in terms of computation cost than 
the coupled cluster method, have been used for much 
larger fullerenes and carbon nanotubes. For example, Ceo 
[25-30] has been studied with self-consistent field (SCF) 
and Moller-Plesset second-order (MP2) theory, C240 [31] 
and carbon nanotubes [32] with density functional the- 
ory (DFT) [33,34], and C540 [35] with the Hartrec-Fock 
(HF) method. 

C60 has the highest symmetry Ih in the point group, 
two types of carbon-carbon (CC) bonds (one single bond 
C-C shared by adjacent five- and six-membered car- 
bon rings and one double bond C=C by adjacent six- 
membered carbon rings), and two kinds of bond angles 



(one angle between two adjacent C-C bonds, and another 
one between a C=C bond and a adjacent C-C bond) 
[7]. In section II, we perform full geometry optimiza- 
tions of C48N12 as well as Cgo with both DFT and re- 
stricted HF (RHF) methods. It is found that the C4 8 Ni 2 
aza-fullerene has several distinguishing features: only one 
nitrogen atom per pentagon, two nitrogen atoms prefer- 
entially sitting in one hexagon, S6 symmetry, six unique 
nitrogen-carbon (NC) and nine unique CC bonds. The 
Mulliken charge analysis shows that the doped nitrogen 
atoms in C 4 sNi 2 exist as electron acceptors and three- 
fourths of the carbon atoms as electron donors. Our 
best CC bond lengths and radius of Cgo calculated with 
B3LYP/6-31G* are in excellent agreement with experi- 
ment. 

Total energy calculations of the optimized C48N12 and 
Ceo are discussed in section III. It is found that the high- 
est occupied molecular orbital (HOMO) is a doubly de- 
generate level with a g symmetry and the lowest unoccu- 
pied molecular orbital (LUMO) is a nondegenerate level 
with a u symmetry. The calculated HOMO-LUMO en- 
ergy gap of C 4 sNi2 is about 1 eV smaller than that of 
Ceo- For both molecules, the total energies calculated 
with STO-3G, 3-21G and 6-31G basis sets differ from 
the 6-31G* results by about 1.5%, 0.6% and 0.05%, re- 
spectively, and the HOMO-LUMO gaps decrease about 5 
cV due to electron correlations. For Ceo, our calculated 
results are in agreement with other groups' calculations, 
and our best HOMO-LUMO energy gap calculated with 
B3LYP/6-31G* is in agreement with experiment. 

When an external electric field is applied to a molecule, 
its charges are redistributed and dipoles are induced or 
reoriented [36] . The relation between the dipole moment 
P and the applied field G can be written as [9] 

P = P + c*G + ^G 2 + |g 3 + ... , (1) 

where G is the electric field, Po is the permanent dipole 
moment, a is the dipole polarizability, (3 is the first-order 
hyperpolarizability, and 7 is the second-order hyperpolar- 
izability. The static dipole polarizability (SDP) measures 
the ability of the valence electrons to find a configuration 
which screens a static external field [9] . It has been shown 
that molecules with many delocalized valence electrons 
should display large SDPs [9,37]. The first- and second- 
order hyperpolarizabilities play a key role in the descrip- 
tion of nonlinear optical phenomena since a time- varying 
polarization can act as the source of new components of 
the electromagnetic field [36] . In section IV, we calculate 
the SDPs and first- and second-order hyperpolarizabili- 
ties of C48N12 and Ceo- In comparison to isolated carbon 
or nitrogen atoms, we find an enhancement in the SDP 
due to the delocalized tt electrons in C 4 sNi2 and Ceo- 
The calculated SDP for C*6o is in agreement with exper- 
iment. The average second-order hyperpolarizability of 
C 4 8Ni2 is about 55% larger than that of Ceo- 

When a material is doped, its mechanical, electronic, 
magnetic and optical properties change [7,9-11]. The 
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ability to control such induced changes is vital to progress 
in material science. Raman and infrared (IR) spectro- 
scopic techniques [38,39] are useful experimental tools to 
investigate how doping modifies the structural and dy- 
namical properties of the pristine material and to under- 
stand the physical origin of such induced changes. Over 
the past 10 years, both techniques have been used widely 
to study the vibrational properties of Cgo [40-45], its 
derivative compounds [46-57], and (doped) carbon nan- 
otubes [58-64]. It has been shown that Ceo has in total 
46 vibrational modes including 4 IR-active [40-43] and 
10 Raman-active vibrational modes [44,45]. These stud- 
ies have offered a good guide to the phonon spectrum 
in the solid state of these materials. In section V, we 
perform a vibrational analysis and calculate the infrared 
(IR) intensities of C48N12. Fifty eight IR-active (i.e., 
29 doubly-degenerate and 29 non-degenerate modes) and 
58 Raman-active (i.e., 29 doubly-degenerate unpolarized 
and 29 non-degenerate polarized) frequencies are deter- 
mined. The best vibrational frequencies and IR results 
for Ceo calculated with B3LYP are in excellent agreement 
with experiment and demonstrate the desirable efficiency 
and accuracy of this theory for obtaining quantitative in- 
formation on the vibrational properties of these materi- 
als. Comparison with other groups' calculations of Ceo 
is made and discussed. 

High resolution NMR [65-67] gives spectra which can 
be analyzed to yield parameters such as the nuclear mag- 
netic shielding a [68] and the nuclear spin-spin coupling 
J [69], which characterize molecular systems and struc- 
tures. Both a and J are determined by the electronic 
environments of the nuclei involved. A satisfactory theo- 
retical description of the distribution of electrons in a 
molecule can lead to reliable predictions of a and J, 
which have a number of applications, such as the identi- 
fication of the conformation or structure of the species 
present in a given sample. In section VI the GIAO 
(gauge-including atomic orbital) [70] and the CSGT (con- 
tinuous set of gauge transformations) [71] methods are 
utilized for calculating the nuclear magnetic shielding 
tensor a in C48N12 and Cgo at both the HF and DFT 
levels of theory. Eight 13 C and two 15 N NMR spectral 
signals are predicted for C48N12. Our best calculated 
NMR results for Ceo are in excellent agreement with ex- 
periment. 

Finally, we end in section VII by giving a summary 
and outlook on the potential applications of C48N12. 



II. OPTIMIZED GEOMETRIC STRUCTURE 

The geometries of both C48N12 and Ceo were fully op- 
timized by using the Gaussian 98 program [72,73], where 
we have employed both RHF and DFT methods. Also we 
discuss the effects of basis sets by considering STO-3G, 
3-21G, 6-31G and 6-31G* [74-79]. For the DFT method, 
we use the B3LYP hybrid functional [80]. 



z 




FIG.l: Geometric structure of C48N12. The site numbers 
{5, 9, 14, 21, 26, 30, 35, 39, 45, 50, 55, 60} are for nitrogen 
atoms, while the others are for carbon atoms. 

In Fig.l, we present the geometry of C48N12. The ab 
initio calculations show that C48N12 has only one nitro- 
gen atom per pentagon and two nitrogen atoms prefer- 
entially sit in one hexagon. The symmetry of C48N 12 
is the Sq point group [19,20]. The optimized distances 
(or radii) Rj from the ith atom to the density center of 
the molecule are listed in Table I. We find that there are 
10 unique radii for C48N12, which suggest that C48N12 
is an ellipsoidal structure and has 10 unique sites (2 for 
N sites and 8 for C sites), while Ceo has an equal ra- 
dius for each carbon atom. As shown in section VI, the 
10 unique sites of C48N12 can be identified by NMR ex- 
periments. Comparing the B3LYP results with the RHF 
results shows that the radii are increased by up to 2% due 
to the electron correlation. Comparing the 6-31G and 6- 
31G*'s results shows that adding polarization functions 
decreases the radius of carbon sites but increases the ra- 
dius of nitrogen sites. In comparison with the results of 
STO-3G and the split valence basis sets, we find that 
increasing the basis size would lead to a decreased ra- 
dius. For Ceo, the radius for each carbon site calculated 
by using B3LYP and the 6-31G* basis set is 3.5502 A, 
which is in excellent agreement with experiment (R = 
3.55 A) [86], the LDA (local density approximation) cal- 
culation (R = 3.537 A) with a pseudopotential approach 
(PPA) [87], and the LDA-based Car-Parrrinello molec- 
ular dynamics (CPMD) simulation (R = 3.55 A ) [88]. 
The success of the B3LYP calculation of the Ceo radius 
demonstrates the importance of electron correlation for 
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an accurate description of a molecular geometrical struc- 
ture. 

The calculated net Mulliken charges Q; of carbon and 
nitrogen atoms in C48N12 are also listed in Table I. There 
are two unique types of nitrogen atoms in the structure. 
The net Mulliken charges Qi/q (q = 1.6 x 10 -19 ) on both 
types of N are negative, for example, —0.5953 C and 
-0.6002 C with B3LYP/6-31G*, -0.7803 C and -0.7829 
C with RHF/6-31G*. The net Mulliken charges of the 
carbon atoms in C48N12 separate into two groups: 1/4 
of carbon atoms with negative Q ; and the remaining 3/4 
with positive Q;. Although the Mulliken charge anal- 
ysis cannot estimate the atomic charges quantitatively, 
their signs can be estimated [89]. From these results, 
we find that the doped nitrogen atoms and one-fourth of 
the carbon atoms exist as electron acceptors, and three- 
fourths of the carbon atoms as electron donors. It should 
be mentioned that we also performed calculations of net 



Mulliken charges of carbon and boron atoms in C48B12 
[90] . We found that the doped boron atoms exist as elec- 
tron donors and all carbon atoms as electron acceptors 
[90]. Therefore, C48N12 and C48B12 have opposite elec- 
tronic polarizations, while Ceo is isotropic. In the case of 
doping into silicon, the V family in the periodic table (for 
example, phosphorous) exists as a donor, while the III 
family (for example, boron) exists as an electron accep- 
tor. Thus, the B- or N-substituted doping in Cgo differs 
greatly from that for silicon. This is due to the unique 
structural and electronic properties of C 60 [7] . With re- 
spect to the electron correlation and the choice of basis 
sets, we find that the absolute value of the net Mulliken 
charge Q; for each atom in C48N12 increases with an in- 
crease of the basis size, but decreases due to the electron 
correlation or by adding polarization functions to a given 
basis set. 



Table I: Net Mulliken charge Qi ( q = 1.6 x 10 19 ) and radius Ri ( 1 A = 0.1 nm ) at the site number n» in C48N12 and Ceo 
calculated by RHF and B3LYP methods with a variety of Pople-style basis sets. 











STO-3G 


3-21G 


6-31G 


6-31G* 


Method 


Fullerene 


Site Number {n;} 


Atom 


Ri 


Qi/9 


R 


Qi/« 


R 


Qi/9 


Ri 


Qi/9 










[A] 


[C] 


[A] 


[C] 


[A] 


[C] 


[A] 


[C] 


RHF 


C48N12 


{ 1, 13, 16, 31, 38, 51} 


C 


3.5452 


0.1073 


3.5059 


0.3856 


3.5082 


0.3546 


3.4988 


0.2926 






{ 2, 12, 29, 32, 37, 52} 


C 


3.5591 


-0.0238 


3.5299 


-0.0897 


3.5271 


-0.0515 


3.5162 


-0.0621 






{ 3, 11, 28, 33, 36, 53} 


C 


3.5602 


-0.0192 


3.5353 


-0.0609 


3.5275 


-0.0344 


3.5389 


-0.0480 






{ 4, 15, 27, 34, 40, 54} 


C 


3.5331 


0.0942 


3.4954 


0.3905 


3.4999 


0.3701 


3.4716 


0.3167 






{ 5, 14, 30, 35, 39, 55} 


N 


3.6974 


-0.2406 


3.5691 


-0.9541 


3.5685 


-0.9552 


3.6151 


-0.7803 






{ 6, 18, 24, 42, 48, 58} 


C 


3.4292 


0.0601 


3.4229 


0.2907 


3.4299 


0.2758 


3.3880 


0.2209 






{ 7, 19, 23, 43, 47, 57} 


c 


3.4023 


0.0893 


3.3874 


0.3363 


3.3983 


0.3058 


3.3409 


0.2962 






{ 8, 20, 22, 44, 46, 56} 


c 


3.4723 


0.0718 


3.4599 


0.3191 


3.4577 


0.3522 


3.4416 


0.2433 






{ 9, 21, 26, 45, 50, 60} 


N 


3.5870 


-0.2344 


3.4985 


-0.9878 


3.4951 


-0.9913 


3.5082 


-0.7829 






{ 10, 17, 25, 41, 49, 59} 


c 


3.4851 


0.0952 


3.4586 


0.3703 


3.4631 


0.3737 


3.4222 


0.3036 




C60 


{1, 2, 3, 4, 5, 6, 60} 


c 


3.5473 





3.5238 





3.5300 





3.5226 





B3LYP 


C48N12 


{ 1, 13, 16, 31, 38, 51} 


c 


3.5871 


0.0827 


3.5302 


0.2770 


3.5365 


0.2147 


3.5171 


0.1961 






{ 2, 12, 29, 32, 37, 52} 


c 


3.5923 


-0.0188 


3.5469 


-0.0422 


3.5458 


-0.0175 


3.5275 


-0.0125 






{ 3, 11, 28, 33, 36, 53} 


c 


3.6019 


-0.0193 


3.5533 


-0.0358 


3.5492 


-0.0173 


3.5395 


-0.0298 






{ 4, 15, 27, 34, 40, 54} 


c 


3.5951 


0.0682 


3.5381 


0.2861 


3.5408 


0.2486 


3.5190 


0.2266 






{ 5, 14, 30, 35, 39, 55} 


N 


3.7175 


-0.1883 


3.5918 


-0.7543 


3.5945 


-0.6690 


3.6187 


-0.5953 






{ 6, 18, 24, 42, 48, 58} 


c 


3.5092 


0.0550 


3.4592 


0.2522 


3.4706 


0.1985 


3.4348 


0.1919 






{ 7, 19, 23, 43, 47, 57} 


c 


3.4882 


0.0679 


3.4328 


0.2635 


3.4465 


0.2103 


3.4049 


0.2124 






{ 8, 20, 22, 44, 46, 56} 


c 


3.5431 


0.0618 


3.4917 


0.2637 


3.4946 


0.2623 


3.4720 


0.1998 






{ 9, 21, 26, 45, 50, 60} 


N 


3.6302 


-0.1814 


3.5336 


-0.7787 


3.5304 


-0.7047 


3.5335 


-0.6002 






{ 10, 17, 25, 41, 49, 59} 


c 


3.5554 


0.0723 


3.5056 


0.2686 


3.5087 


0.2548 


3.4807 


0.2110 




C60 


{1, 2, 3, 4, 5, 6, 60} 


c 


3.6034 





3.5555 





3.5615 





3.5502 
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Table II: Bond lengths (L, 1 A = 0.1 nm ) in C48N12 and Ceo calculated by using B3LYP methods with a variety of Pople-style 
basis sets, where (m,nj) denotes the site number pair that forms a bond. 







STO-3G 


3-21G 


6-31G 


6-31G* 


FuUerene Bond 




Ld/t 




Ld/t 


L r h/ 


Ld/t 


Lrh/ 


Ld/t 


L r h/ 






[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


[A] 


C48N12 CC 


(1 , 2) (12, 13) (16, 29) 
(31, 32) (37, 38) (51, 52) 


1.4275 


1.3914 


1.4103 


1.3855 


1.4125 


1.3880 


1.4061 


1.3836 


CC 


(1 , 52) (2 , 31) (12, 38) 
(13, 29) (16, 37) (32, 51) 


1.4371 


1.4151 


1.4171 


1.4083 


1.4216 


1.4134 


1.4155 


1.4024 


NC 


(1 , 5) (13, 14) (16, 30) 
(31, 35) (38, 39) (51, 55) 


1.4742 


1.4601 


1.432 


1.4206 


1.4315 


1.4164 


1.4300 


1.4272 


CC 


(2, 3) (11, 12) (28, 29) 
(32, 33) (36, 37) (52, 53) 


1.4701 


1.4662 


1.4517 


1.4524 


1.4488 


1.4478 


1.4455 


1.4521 


CC 


(3, 4) (11, 15) (27, 28) 
(33, 34) (36, 40) (53, 54) 


1.4092 


1.3590 


1.3947 


1.3620 


1.3971 


1.3652 


1.3901 


1.3588 


CC 


(3, 10) (11, 59) (17, 33) 
(25, 36) (28, 49) (41, 53) 


1.4559 


1.4566 


1.4320 


1.4350 


1.4346 


1.4367 


1.4314 


1.4368 


NC 


(4 , 5) (14, 15) (27, 30) 
(39, 40) (34, 35) (54, 55) 


1.4707 


1.4553 


1.4307 


1.4128 


1.4286 


1.4108 


1.4224 


1.4047 


CC 


(4, 46) ( 8, 15) (20, 40) 
(22, 54) (27, 44) (34, 56) 


1.4567 


1.4589 


1.4326 


1.4367 


1.4354 


1.4388 


1.4313 


1.4372 


NC 


(5, 42) ( 6, 35) (14, 48) 
(18, 55) (24, 30) (39, 58) 


1.4749 


1.4637 


1.4299 


1.4294 


1.4317 


1.4278 


1.4287 


1.4310 


CC 


(6 , 7) (18, 19) (23, 24) 
(42, 43) (47, 48) (57, 58) 


1.4399 


1.4352 


1.4211 


1.4182 


1.4217 


1.4185 


1.4136 


1.4143 


CC 


(6 , 10) (17, 18) (24, 25) 
(58, 59) (41, 42) (48, 49) 


1.4119 


1.3641 


1.3963 


1.3664 


1.4004 


1.3712 


1.3941 


1.3634 


CC 


(7,8 ) (19, 20) (22, 23) 
(43, 44) (46, 47) (56, 57) 


1.4195 


1.3705 


1.4021 


1.3718 


1.4068 


1.3761 


1.4021 


1.3703 


NC 


(7, 60) (9, 47) (19, 26) 
(21, 57) (23, 45) (43, 50) 


1.4586 


1.4419 


1.4159 


1.4047 


1.4189 


1.4071 


1.4099 


1.3958 


NC 


(8 , 9) (20, 21) (22, 26) 
(44, 45) (46, 50) (56, 60) 


1.4510 


1.4313 


1.4182 


1.4083 


1.4149 


1.4012 


1.4084 


1.4056 


NC 


(9 , 10) (17, 21) (25, 26) 
(41, 45) (49, 50) (59, 60) 


1.4612 


1.4360 


1.4271 


1.4043 


1.4226 


1.3997 


1.4134 


1.3919 


Ceo c=c 


(1, 52) (2, 31) (3, 10) ... 


1.4130 


1.3759 


1.3899 


1.3671 


1.3981 


1.3750 


1.3949 


1.3732 


C-C 


(1, 2) (1, 5) (2, 3) ... 


1.4773 


1.4628 


1.4601 


1.4529 


1.4592 


1.4524 


1.4539 


1.4487 



The optimized CC and NC bond lengths in C48N12 are 
listed in Table II. We find that there are 6 unique NC and 
9 unique CC bonds. In comparison with the calculated 
CC bond lengths for Ceo shown in Table II, the CC bond 
length in C4gN 12 , in general, is less than the single C- 
C bond length of Cgo due to the redistribution of the 
electron density. It is also found that the bond length 
increases due to the electron correlation, but decreases 
as we increase the basis size or include the polarization 
function. 

In comparison with experimental data available for Ceo 
listed in Table III, we find that the two kinds of bond 
lengths of C@o calculated by using the B3LYP with a 
large basis set 6-31G* are in good agreement with the re- 
sults measured by X-ray powder diffraction (XRPD) [91], 
NMR [92,93], gas-phase electron diffraction (GPED) [94] 
or X-ray crystallography technique (XRCT) [86]. 



For comparison, Table III also lists the calculated 
CC bond lengths for C@o with a selection of previous 
theoretical calculations. Given the low computational 
cost of Hiickcl theory, the bond lengths [95] predicted 
by this theory are remarkably satisfactory. The semi- 
empirical QCFF/PI (quantum-chemical- force- fields for it 
electrons) [96] does not predict as good bond lengths as 
the Hiickel theory since it has been parameterized mainly 
with respect to frequencies of conjugated hydrocarbons. 
The CC bond lengths calculated by using the semi- 
empirical MNDO (modified neglect of differential over- 
lap) [97] and the extended Hubbard model (EHM) [98] 
are a little improved. These theoretical approaches em- 
pirically include the effect of electron correlation found in 
conjugated 7r-systems. The HF [25-27,30] and SCFMO 
(self-consistent field with MO) [29] calculations are in 
agreement with our RHF results. As listed in Table III, 
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the calculated MP2 bond distances [28] usually decrease 
by about 0.01 A when more d functions are added to 
the basis set, demonstrating the necessity of including 
polarization functions in calculations with correlation. 
Based on the differences between the HF and MP2 data, 
it is evident that electron correlation effects should be 
considered in an accurate description of the equilibrium 
structure of a molecule. The CC bond lengths calcu- 
lated with the LDA [99-101], LDA-PPA [87], and LDA- 
based CPMD simulation [88,102] are in agreement with 
our B3LYP's results and demonstrate the importance of 
electron correlation effects in giving accurately the equi- 
librium structure of a molecule. 

Table III: Equilibrium CC bond lengths (1 A = 0.1 nm) 
of C6o from previous theoretical predictions and experimen- 
tal findings. STO-DZP and DNP denote double-zeta STO's 
and double numerical basis with polarization functions, re- 
spectively. 



Method 


C-C 




Reference 




[A] 


[A] 




Huckel 


1.436 


1.418 


[951 


QCFF/PI 


1.471 


1.411 


[96] 


MNDO 


1.465 


1.376 


[97] 


EHM 


1.446 


1.402 


[98] 


SCFMO 


1.49 


1.43 


[29] 


HF/STO-3G 


1.463 


1.376 


[25] 


HF/3-21G 


1.453 


1.367 


[30] 


HF/DZ 


1.451 


1.368 


[26] 


HF/STO-3G 


1.463 


1.376 


[27] 


HF/DZ 


1.451 


1.368 


[27] 


HF/DZP 


1.450 


1.375 


[27] 


HF/TZP 


1.448 


1.370 


[27] 


MP2/DZ 


1.470 


1.407 


[28] 


MP2/DZP 


1.451 


1.412 


[28] 


MP2/TZP 


1.446 


1.406 


[28] 


LDA/STO-DZP 


1.436 


1.384 


[101] 


LDA/DZP 


1.445 


1.395 


[100] 


LDA/DNP 


1.444 


1.391 


[99] 


LDA-PPA 


1.449 


1.390 


[87] 
[88] 


LDA-CPMD 


1.45 


1.40 


LDA-CPMD 


1.45 


1.39 


[102] 


Exp./XRCT 


1.4459 


1.3997 


[86] 


Exp./XRPD 


1.455 


1.391 


[91] 


Exp./NMR 


1.45 


1.40 


[92] 


Exp./NMR 


1.46 


1.40 


[93] 


Exp./GPED 


1.458 


1.401 


[94] 



As mentioned before, Ceo has only two kinds of bond 
angles [7], 108° (the angle between two adjacent sin- 
gle C-C bonds) and 120° (the angle between a double 
C=C bond and an adjacent single C-C bond). Fig.2(a-d) 



show the distribution of (C-C-C, C-N-C, C-C-N) bond 
angles in C4gN 12 calculated by using both RHF and 
B3LYP methods with several different basis sets. One- 
third and two-thirds of the bond angles in C48N12 fluc- 
tuate around 108° and 120°, respectively. Comparing 
Fig. 2 (a) and 2(b) shows that increasing the basis size 
leads to smaller fluctuations in the bond angle distribu- 
tion (BAD) around either 108° or 120°. In comparison 
with Fig. 2(b), Fig. 2(c) exhibits enhanced fluctuations in 
the BAD as polarization functions are added to the 6- 
31G basis set. Comparing Fig. 2(c) with Fig. 2(d), we find 
that the BAD in C48N12 are decreased due the effect of 
electron correlation. 



III. TOTAL ELECTRONIC ENERGY 

We performed total energy calculations of C48N12 and 
C 60 by using both RHF and B3LYP with STO-3G, 3- 
21G, 6-31G and 6-31G* basis sets. The results are sum- 
marized in Table IV. The orbital energies are shown in 
Fig. 3, where the orbital symmetries are also labeled. 

Table IV demonstrates the convergence of the total 
energy calculations of both RHF and B3LYP methods 
with respect to the basis sets. For both C48N12 and Ceo, 
the total energies calculated with STO-3G, 3-21G and 
6-31G basis sets differ from the 6-31G* basis set results 
by about 1.5%, 0.6% and 0.05%, respectively. For both 
molecules, comparing the B3LYP and RHF results shows 
that the HOMO-LUMO energy gap A decreases about 5 
eV and the binding energy Eb per atom increases about 
2 eV because of the electron correlations. The calculated 
binding energy Eb per atom and HOMO-LUMO energy 
gap A of C48N12 are about 1 eV smaller than those of 

C60- 

Because of the valency of the doped nitrogen atoms, 
the electronic properties of C48N12 and Ceo are signifi- 
cantly different. As shown in Fig. 3(a) and Fig. 3(b), the 
HOMO for Ceo is fivefold-degenerate with h u symmetry, 
the LUMO is threefold-degenerate with ti u symmetry, 
and the others are threefold-degenerate with t\ g or t2 U 
symmetry, fourfold-degenerate with g g or g u symmetry, 
and fivefold-degenerate with h g symmetry. We notice 
from Fig. 3 that each energy level in C48N12 splits since 
the icosahedral symmetry of Ceo is lost by the substi- 
tutional doping. For C48N12, the HOMO is a doubly 
degenerate level of a g symmetry, the LUMO is a non- 
degenerate level with a u symmetry, and the others are 
specified in Fig. 3. Considering that C48N12 is isoelec- 
tronic with Cg 12 , we find that the filling of the energy 
levels in C48N12 corresponds to a complete filling of the 
ti u and tig levels of Cgo- 
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FIG. 2: Ab initio calculations of C-C-C (open circles) and C-N-C (or C-C-N) (filled circles) bond angles in C48N12: 
(a) B3LYP/STO-3G; (b) B3LYP/6-3fG; (c) B3LYP/6-3f G*; (d) RHF/6-31G*. The site numbers are labeled as 
in Fig.f. The solid lines are the C-C-C bond angles ( f08° and 120°) of C 6 o- 



Table IV: Total electronic energy (E t , in eV), LUMO energy (E lumo , in eV), HOMO energy (E homo , in eV), HOMO-LUMO 
energy gap (A, in eV), binding energy per atom (Et,, in eV), ionization potential (Eip, in eV) and electron affinity (Eea, in 
eV) of C48N12 and C60 calculated by using RHF and B3LYP methods with a variety of Pople-style basis sets. 



Method 


Energy 




C60 








C 48 N 


12 




STO-3G 


3-21G 


6-31G 


6-31G* 


STO-3G 


3-21G 


6-31G 


6-31G* 


RHF 


E t 


-61067.132 


-61471.413 


-61795.062 


-61817.394 


-66397.387 


-66846.232 


-67193.956 


-67223.547 






2.688 


-0.632 


-0.709 


-0.117 


3.797 


0.499 


0.214 


0.643 




Ehomo 


-5.456 


-8.323 


-7.919 


-7.644 


-3.287 


-6.506 


-6.203 


-6.189 




A 


8.144 


7.691 


7.210 


7.527 


7.084 


7.005 


6.417 


6.832 




E b 


5.581 


4.627 


4.666 


4.956 


4.510 


3.734 


3.724 


4.149 




Eip 


4.892 


8.230 


5.685 


5.659 


2.45 


6.144 


5.314 


5.530 




Eea 


2.280 


0.881 


1.187 


1.125 


1.947 


0.146 


0.209 


0.267 


B3LYP 


E t 


-61446.662 


-61864.864 


-62194.006 


-62209.002 


-66795.835 


-67263.771 


-67617.031 


-67637.724 




Eluino 


-1.127 


-3.563 


-3.390 


-3.219 


-0.188 


-2.783 


-2.706 


-2.614 




Ehomo 


-4.356 


-6.509 


-6.221 


-5.987 


-2.365 


-4.743 


-4.633 


-4.383 




A 


3.229 


2.946 


2.831 


2.768 


2.177 


1.96 


1.927 


1.774 




E b 


7.715 


6.786 


6.843 


6.982 


6.888 


6.140 


6.088 


6.365 




Eip 


3.868 


7.814 


8.576 


7.317 


2.488 


6.040 


5.600 


5.663 




Eea 


1.487 


2.295 


1.073 


2.398 


0.804 


1.514 


0.223 


1.493 
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In Table IV, we list the calculated ionization poten- 
tial (IP) Eip and electron affinity (EA) Eea for both 
Ceo and C48N12. It shows that C48N12 is a good elec- 
tron donor, while Cgo is a good electron acceptor. The 
calculated IP for Cgo is in good agreement with the ex- 
periments: (7.54 ±0.01) eV [103], (7.57 ±0.01) eV [104], 
(7.58±g;o2 eV [105], and (7.59 ± 0.02) eV [106]. The cal- 
culated EA for Ceo agrees well with the experiments: 
(2.666 ±0.001) eV [107], (2.689 ± 0.008) eV [108]. 
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Fig. 3: Orbital energies Efe of the fcth eigenstate 
of Ceo and C48N12 calculated with (a) B3LYP/6- 
31G* and (b) RHF/6-31G*. Open circles are for 
Ceo- Filled circles and triangles are for C48N12, 
respectively. The orbital symmetries of energy 
levels are shown. 

For comparison, here we list the total energies of Ceo 
for a selection of previous theoretical calculations. Us- 
ing the HF methods, Scuseria [27] found that the total 
energies of Ceo calculated with DZ, DZP and TZP ba- 
sis sets [25-27,30] are -61745.723 eV, -61826.192 eV and 
-61832.443 eV, respectively. These results are in agree- 
ment with our RHF calculations. Moreover, using the 
MP2 method, Scuseria and colleagues [28] have shown 
that the total energies of Ceo with DZ, DZP and TZP ba- 
sis sets are -61893.398 eV, -61826.192 eV and -61832.443 
eV, respectively. The total energies calculated by using 
the MP2 method do improve the approximation made 
in HF method. This result demonstrates that electron 
correlation effects cannot be neglected for an accurate 
prediction of the total energy of a molecule. 



In crystals, the on-site Coulomb interaction between 
two electrons on the same molecule is given by [109] 



U — Ufree — 2E p , 
Ufree = E/p — EeA — A, 



(2) 
(3) 



where E p = zae 2 / (2L 4 ) (z: the number of nearest neigh- 
bors, L: the distance between molecules, a: the dipole 
polarizability) is the polarization energy and Uf ree is 
the on-site Coulomb interaction between two electrons 
on the free molecule [7]. Our plane- wave pseudopoten- 
tial calculations performed by using the CASTEP pro- 
gram [73,110] show that the optimized lattice constants 
for Cgo and C4sNi2-based fee solids (note: z = 12 for 
both cases) are all around 1.45 nm, while the experimen- 
tal lattice constant = 1.4161 nm for Ceo [HI]- Based on 
our first principles results for Eip, Eea and A for Ceo 
and C48N12 molecules, we arrive at the value of the on- 
site Coulomb interaction, U = 1.84,2.05 eV for Ceo and 
C48N12 solids, respectively. The experimental values for 
C 60 are U = 1.6 ± 0.2 eV [109] and 1.54 eV [112]. More- 
over, our DFT/GGA calculations show that the band gap 
Egap for the C48N12 solid is about 0.85 eV smaller than 
that of Ceo- This is consistent with the A-U-W relation 
[7] 



E, 



gap 



A + U-W, 



(4) 



where W is the bandwidth for the HOMO- or LUMO- 
derived energy bands. Assuming the same W for both 
C48N12 and Cgo solids, we arrive at 



C 48 Ni 



gap 



z?C eo 

gap 



0.79 eV. 



(5) 



The DFT is known to underestimate the band gap of 
solids. From experiment, E gap = 2.3 ±0.1 eV [109] and 
2.86 eV [112] for the Ceo solid. Hence, using the approx- 
imate relation given by Eq.(5), we obtain an estimated 
band gap for the C 48 Ni2 fee solid of E^ Nl2 = 1.7 ± 0.3 
eV. Thus, the C48N12 solid, like Ceo, is a semiconducting 
material. 

Similarly, we find that the C48Bi 2 -based fee solid (lat- 
tice constant w 1.45 nm) is also a semiconducting ma- 
terial [90], having U = 1.9 eV, E^f^ - = -1.34 
eV, and E^ Bl2 = 1.2 ± 0.3 eV. 



IV. STATIC (HYPER)POLARIZ ABILITY 

The static polarizabilities for C48N12 and Ceo are pre- 
sented in Table V. The B3LYP and RHF results are ob- 
tained by using the Gaussian 98 program package [72,73] 
and the LDA results by using the ADF (Amsterdam Den- 
sity Functional) program [73,113,114]. The ADF pro- 
gram uses basis sets of Slater functions, where a triple 
zeta valence basis plus polarization is augmented with 
the field-induced polarization (FIP) functions of Zeiss 
ct al. [115]. Here this basis set is denoted as TZP±± 
([6s4p2dlf] for C and N atoms, and [4s2pld] for H atom) 
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and has recently been used for calculating the second- 
order hyperpolarizabilities, 7, for Ceo, C58N2, C58B2 and 
C 58 BN [116]. 

From the results listed in Table V, we see that the ba- 
sis set dependence is identical for both B3LYP and RHF 
cases. As expected, improving the basis set increases the 
polarizability. Our RHF results for Ceo are in good agree- 
ment with a previous study [117] which used a similar 
method and basis sets. The B3LYP values are about 10% 
larger than the corresponding RHF values. The LDA re- 
sults are much larger than both the B3LYP and RHF 
results. This is expected since the basis set is larger 
and expected to predict a more accurate polarizability 



[116]. Also LDA, in general, predicts a larger polariz- 
ability than docs with B3LYP [118]. The LDA results 
for Ceo are in good agreement with previous LDA stud- 
ies [119-121]. We find that the polarizability of C48N12 
is slightly smaller than the polarizability of Ceo- For 
C48N12 we also find that the zz component is slightly 
larger than the xx component except when using RHF 
and a small basis set. In comparison with the SDP of 
a single carbon or nitrogen atom, we find an enhanced 
linear polarizability for both Ceo and C48N12. This is 
expected since both molecules have many conjugate tt 
electrons dclocalized over the entire system. 



Table V: Static polarizabilities (a, in A 3 = 10 30 m 3 ) for C4SN12 and C60 with RHF and DFT methods and a variety of 
basis sets. Symmetry relations give for C and for C48N12, O-xx — ^yy~ 





STO-3G 




3-21G 




6-31G 




6-31G* 




TZP++ 




B3LYP 


RHF 


B3LYP 


RHF 


B3LYP 


RHF 


B3LYP 


RHF 


LDA 


Ceo 

Oi-xx 


52.5 


45.2 


64.9 


59.5 


68.3 


63.2 


69.5 


64.7 


84.7 


C48N12 

&XX 


51.5 
51.6 


43.5 
41.6 


62.0 
62.6 


55.6 
54.9 


65.6 
66.1 


59.4 
58.5 


66.6 
67.5 


60.0 
60.4 


79.3 
81.5 



Table VI: Static 2nd-order hyperpolarizabilities (7, in au. 1 au = 6.235378 xlCT 65 C 4 m 4 J- 3 ) for C48N12 and C 60 
calculated using LDA and TZP++. The average 2nd-order hyperpolarizability is given by 7 = J2i j (lajj + Jijij + 
■Jijji). Symmetry relations give j xxxx = j yy yy, J xxzz = J vyzz and j zzxx = j zzyv . 





"fxxxx 


fxxyy 


'Jzzzz 


'Jxxzz 


^fzzxx 


7 




137950 


45983 


137950 


45983 


45983 


137950 


C48N12 


188780 


62880 


232970 


85120 


84790 


215222 



Because Cqo is well separated in a crystal [7], one can 
model the fullcrcnc crystal in an electric field as a collec- 
tion of isolated dipoles. For such systems, the Clausius- 
Mossotti relation [122] is expected to yield a reasonably 
accurate relationship between the linear polarizability 
of an isolated molecule and the dielectric constant of a 
fullerene crystal, i.e., 



a = 



3(e~l) 
47rp/(e + 2)' 



(6) 



where pf is the density of fullerene molecules in a fee 
crystal and e is the dielectric constant of the fullerene 
crystal. For the purposes of comparison, we note that 
Hebard et al. [123] and Ecklund [124] measured the di- 
electric constants of Cgo films in the range between 3.9 
and 4.0. Based on these measured dielectric constants, 
one can derive an experimental polarizability of about 
84.9 x 10~ 30 m 3 for C 60 (the density p f = 10 30 /712 m~ 3 
is given by Quong and Pederson [125]). Very recently, us- 
ing the molecular beam deflection technique, Antoine et 



al. [126] have measured the electric polarizability of iso- 
lated Cgo molecules and obtained a value of (76.5 ±8.0) x 
10~ 30 m 3 . Using a new optical technique that uses light 
forces and a time-of-fiight spectrometer, Ballard et al. 
[127] have made absolute measurements of cluster po- 
larizabilities and determined the optical polarizability of 
Ceo at the fundamental wavelength of a NdYAG laser 
( A = 1064 nm) to be (79 ± 4) x 10~ 30 m 3 . The ex- 
perimental results are in good agreement with our LDA 
results, especially considering that LDA is expected to 
overestimate the polarizability. 

The static first-order hyperpolarizability (3 [36] of the 
C48N12 molecule is also calculated and found to be zero, 
the same as that for Ceo molecule. This is expected since 
both C48N12 and Gqq molecule display inversion symme- 
tries. Consequently, this aza-fullerene cannot produce 
second-order nonlinear optical interactions. 

The static 2nd-order hyperpolarizabilities for C48N12 
and Cgo molecules are presented in Table VI. For 
the calculations of the 2nd-order hyperpolarizability, 
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7, we use time-dependent (TD) DFT as described in 
Rcf. [116,120,128,129]. First, the lst-order hyperpolar- 
izability, (3, is calculated analytically in the presence of 
a small electric field. Then, the 2nd-ordcr hypcrpolar- 
izability can be obtained by a finite-field differentiation 
of the analytically calculated lst-order hyperpolarizabil- 
ity. For all the TD-DFT calculations we used the RE- 
SPONSE code [73,128,130] implemented in the ADF pro- 
gram [73,113,114]. The small difference between "f xxzz 
and jzzxx for C48N12 is due to the numerical method 
adopted for calculating 7. 

For the 7 value of Ceo, we find good agreement with 
previous first-principles results [120,121,131]. A compar- 
ison with experiments will not be made for the 2nd-order 
hyperpolarizability due to large differences in the exper- 
imental results [9,10,120,131]. We find that all compo- 
nents of the 2nd-order hyperpolarizability for C48N12 
are larger than for Cw This gives an average 2nd- 
order hyperpolarizability of C48N12 which is about 55 
% larger than the average 2nd-order hyperpolarizability 
of Ceo- The zzzz components of the 2nd-order hyperpo- 
larizaility of C48N12 is also larger than that found in the 
donor/acceptor substituted C58BN molecule [116]. 



V. VIBRATIONAL FREQUENCY ANALYSIS 
A. Theory 

To lowest order, IR intensities are proportional to the 
derivatives of the dipolc moment with respect to the vi- 
brational normal modes of the material, evaluated at the 
equilibrium geometry. In detail, the IR intensity of the 
qth vibrational mode is given by [132] 



dPj 

dR k 



d 2 E 
'dGidR k 



(9) 



(q) _ PpK 

IR ~ 3c 



dP 



(7) 



where p p is the particle density, S g is the normal coordi- 
nate cooresponding to the qth mode and c is the velocity 
of light. Since \dP/dE q \ is the only molecular property 
entering the formula, it is often referred to as absolute 
IR intensity. 

To obtain the IR data, one must compute the deriva- 
tives of the dipole moment with respect to the normal 
mode coordinates. These can be viewed as directional 
derivatives in the space of 3N nuclear coordinates and 
expressed in terms of derivatives with respect to atomic 
coordinates, Rk- For the ith component of the dipole 
moment P (i=x,y,z), we have 



dPj 

dZ„ 



3N dPi 



fe=i 



dRk 



(8) 



where £,k q — dRk/dE q is the fcth atomic displacement of 
the gth normal mode. Then, the necessary derivatives 
can be expressed in terms of the atomic forces as follows 
[132,133] 



where E is the total energy, d is the ith component of 
an assumed external electric field G, and Fk is the cal- 
culated force on the fcth atomic coordinate. 



B. Normal Vibrations in C60 and C48N12 

Using the Gaussian 98 program [72,73], we calculated 
the harmonic vibrational frequencies of both Cgo and 
C48N12 and considered the effects of the basis sets. It 
should be mentioned that our frequencies have not been 
scaled. 

Table VII and VIII summarize the vibrational frequen- 
cies for Ceo calculated by using RHF and B3LYP meth- 
ods, respectively. As shown by Dresselhaus et al. [134], 
there are 46 vibrational modes for Cgo- These modes are 
classified in even and odd parities and in ten irreducible 
representatives of the Ih point group: the {a g , a„}, {ti g , 
tin, h g , t 2u }, {g g , 9u} and {h g , h u } modes are non-, 
threefold-, fourfold- and fivefold-degenerate, respectively. 
Tables VII and VIII demonstrate that increasing the ba- 
sis size improves the accuracy of the predicted vibrational 
frequencies, but adding polarization functions to the 6- 
31G basis set only improves slightly the accuracy of the 
vibrational frequencies. In choosing a basis set for the 
first-principles calculation, one must make a compromise 
between accuracy and CPU time. Our results shows that 
the minimum calculation can be done in about 4 hours 
of CPU, while the most expensive calculation requires 
12 days of CPU. Without significant computational cost, 
one can do B3LYP/STO-3G calculation and still obtain 
results more accurate than any RHF calculations. Go- 
ing beyond STO-3G for B3LYP calculations requires a 
drastic increase in CPU time. Suprisingly, going just to 
3-21G provides the most accurate results, while for the 
bigger basis set 6-31 G, the results are worse and adding 
a polarized function to 6-3 1G only slightly improves the 
results. The 3-21G basis set gives systematically lower 
frequencies than the 6-31G basis set, while the frequen- 
cies obtained from the 6-31G* basis set typically lie be- 
tween the results of the other two basis sets. In contrast, 
as discussed in section II, 6-31G* does provide the most 
accurate bond lengths. This suggests that the better ac- 
curacy of 3-21G is foirtuitous. Increasing the basis set to 
6-31G stiffens the bonds, while adding the polarization 
function compensates by softening the bonds. In compar- 
ison with the B3LYP results, RHF calculated frequencies 
are too high due to an incorrect description of bond dis- 
sociation, while B3LYP with large basis sets (even the 
minimum basis set STO-3G) generally gives results in 
good agreement with the experiments of Wang et al. [42] 
and Dong et al. [47]. This demonstrates the importance 
of electron correlation in an accurate description of the 
vibrational frequencies. 
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For comparison, Table VII lists the vibrational fre- 
quencies of Ceo calculated by using various theories, for 
example, the semi-empirical MNDO [97] and QCFF/PI 
[96] methods. Of these, the QCFF/PI method, which 
has been parameterized mainly with respect to vibra- 
tional frequencies of conjugated and aromatic hydrocar- 
bons [135], results in the best results although it gives less 
satisfactory geometry. Such accurate prediction implies 
that the electronic structures of Ceo is not much differ- 
ent from other aromatic hydrocarbons [96]. Haser et al. 
[28] showed that the approximate harmonic frequencies 
for the two a g vibrational modes of Cgo are 1615 cm -1 
and 487 cm" 1 for HF/DZP, 1614 cm" 1 and 483 cm" 1 for 
HF/TZP, 1614 cm" 1 and 437 cm" 1 for MP2/DZP, and 
1586 cm" 1 and 437 cm" 1 for MP2/TZP. Their HF cal- 
culations are in agreement with our RHF/3-21G results. 
Their MP2 results are more accurate when obtained with 
large basis sets, which also demonstrates the importance 
of electron correlation in predicting the vibrational fre- 
quencies. 

In addition, there have been a number of 2nd nearest- 
neighbor force-constant models (FCMs) [136-138] which 
have been used to calculate the phonon frequencies of 
Ceo- None of them yield good agreement with the ex- 
perimental data. For example, an empirical force field, 
which has been parameterized with respect to polycyclic 
aromatic hydrocarbons, is used with Hiickel theory and 
predicts vibrational frequencies of the two a g modes of 
1409 cm -1 and 388 cm" 1 [137] that are too low. However, 
the modified FCM (MFCM) by Jishi et al. [139] consid- 
ered interactions up to the third-nearest neighbors, and 
the calculated results, as shown in Table VII, are in ex- 
cellent agreement with the experiments of Wang et al. 
[42] and Dong et al. [47]. 

Table VIII also lists the vibrational frequencies of Ceo 
calculated by other DFT methods, for example, LDA- 
PPA [87], LDA [100,101] and DFT-LDA-based CPMD 
simulations [102]. In general, those calculated results are 
in good agreement with experiment. Very recently, Choi 
et al. [141] have performed B3LYP vibrational calcula- 
tions of Ceo with a 3-2 1G basis set but involving scaling 
of the internal force constants (SIFC) K™ 1 by using Pu- 
lay's method [142], i.e., 

^scaled = (s . s ^l/2^mt ; (1Q) 

where iC™* is the force constant in internal coordinates 
( the Gaussian 98 program [72] uses this form), and Sj 
and Sj are scaling factors for the ith and jth redundant 
internal coordinates, respectively. They optimized the 
scaling factors by minimizing the root-mean-square de- 
viations between the experimental and calculated scaled 
frequencies. Their results are listed in Table IX. Overall, 
their scaling procedure improves the accuracy for the 46 
vibrational frequencies of Ceo, especially, for the a g , h g 
and t\ u vibrational modes. 

In Table X and XI, we list the vibrational frequencies 
for C 48 Ni2 calculated with RHF and B3LYP methods 



and a variety of Pople-style basis sets. In contrast with 
Ceo, it is found that there are in total 116 vibrational 
modes for C48N12 because of its lowered symmetry, Sq. 
These vibrational modes are classified into 58 doubly- 
degenerate and 58 nondegenerate modes. Among those 
vibrational modes, there are 58 IR-active (listed in Ta- 
ble X) and 58 Raman-active modes (listed in Table XI). 
Table X and XI show that the electron correlation or 
increasing the basis size results in a redshift of the vibra- 
tional frequencies. This is similar to that of Ceo- 



C. IR Intensities in C48N12 and C60 

We perform calculations of IR intensities Ijh for both 
C48N12 and Ceo by using the Gaussian 98 program [72,73] 
with RHF and B3LYP methods. The calculated IR inten- 
sities for Ceo at the corresponding frequencies are listed 
in Table XII, and those for C48N12 are shown in Fig. 4. 

For Ceo, we note that its IR spectrum is very simple. 
Namely, it is composed of 4 IR-active vibrational modes 
with ti u symmetry. This is a consequence of the sym- 
metry of the icosahcdral group [134]. Carbon clusters of 
comparable size, but lower symmetry, have many more 
IR-active frequencies. For example, the graphitene iso- 
mer of Ceo has D 6 h symmetry and 20 IR-active frequen- 
cies [97]. Other examples [97] included C54, a planar 
graphite fragment with D 6 h symmetry and 22 IR-active 
frequencies, and C50, a spheroidal cluster with D^t sym- 
metry and 22 IR-active frequencies. From Table XII, 
we see that the IR intensity of a given mode decreases 
due to the electron correlation and converges with in- 
creasing basis size. We find that our intensities calcu- 
lated with B3LYP agree reasonably with experimental 
spectrum [140] obtained by in situ high-resolution FTIR 
measurement of a Ceo film. 

The IR intensities Iir at the corresponding vibrational 
frequencies for C48N12 are presented in Fig.4(a)(b)(c) for 
calculations with B3LYP/STO-3G, B3LYP/3-21G and 
RHF/3-21G, respectively. As discussed above, C48N12 
has 29 nondegenerate and 29 doubly-degenerate IR- 
active vibrational modes. Similar to the case of C@o, the 
IR signals separate into two regions, i.e., a high-frequency 
( > 1000 cm -1 ) and a low- frequency ( < 1000 cm -1 ) 
region. The IR-active frequencies are redshifted by in- 
cluding the electron correlation or increasing the basis 
size. The IR intensities decrease after including elec- 
tron correlations and converge with increasing basis size. 
The strongest IR spectral lines in both low- and high- 
frequency regions are the doubly-degenerate modes lo- 
cated, for example, at 440 cm -1 and 1310 cm -1 , respec- 
tively, for the B3LYP/3-21G case. Since experimental 
IR spectroscopic data do not directly indicate the spe- 
cific type of nuclear motion producing each IR peak, we 
do not give here the normal mode information for each 
vibrational frequency and the displacements of each nu- 
clei corresponding to each normal mode. In Fig. 5, taking 
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B3LYP/3-21G calculations as an example, we only show and high-frequency regions. It is seen that the pentagon 
the vibrational displacements of sites 1 to 5 (4 C sites and structure for the low- frequency case expands slightly and 
1 N site) for the strongest IR spectral signals in the low- 
Table VII: Vibrational frequencies (v, in cm -1 ) of Cgo calculated with RHF and a variety of Pople-style basis sets. Numbers 

in the parenthesis are the relative errors to the experimental frequencies, v exp , from Wang et al. [42] and Dong et al. [47]. 
The approximated CPU times for STO-3G, 3-21G, 6-31G and 6-31G* basis sets are 5 hours, 12 hours, 18 hours and 39 hours, 

respectively. Results of other theoretical calculations, for example, QCFF/PI [96], MNDO [97] and MFCM [139], are also listed. 

Mode STO-3G 3-21G 6-31G 6-31G* QCFF/PI MFCM MNDO Exp = 
Even Parity 

cig 1684 (14.5%) 1604 (9.1%) 1637 (11.3%) 1600 (8.9%) 1442 (1.9%) 1468 (0.1%) 1667 (13.4%) 1470 

553 (11.0%) 518 (4.0%) 526 (5.7%) 527 (5.7%) 513 (3.0%) 492 (1.2%) 610 (22.5%) 498 

g g 1802 (18.2%) 1667 (9.3%) 1697 (11.3%) 1687 (10.6%) 1585 (3.9%) 1521 (0.3%) 1650 (8.2%) 1525 

1531 (12.9%) 1426 (5.1%) 1450 (6.9%) 1441 (6.3%) 1450 (6.9%) 1375 (1.4%) 1404 (3.5%) 1356 

1203 (11.8%) 1095 (1.8%) 1142 (6.2%) 1132 (5.2%) 1158 (7.6%) 1056 (1.9%) 1235 (14.8%) 1076 

908 (12.6%) 787 (2.3%) 878 (9.0%) 836 (3.8%) 770 (4.5%) 805 (0.1%) 856 (6.2%) 806 

653 (5.2%) 646 (4.1%) 643 (3.5%) 619 (0.3%) 614 (1.1%) 626 (0.8%) 579 (6.9%) 621 

574 (18.0%) 529 (9.0%) 549 (12.9%) 537 (10.5%) 476 (2.1%) 498 (2.5%) 491 (1.0%) 486 

h g 1912 (21.2%) 1772 (12.3%) 1799 (14.0%) 1791 (13.5%) 1644 (4.2%) 1575 (0.2%) 1722 (9.1%) 1578 

1658 (16.2%) 1546 (8.4%) 1585 (11.1%) 1562 (9.4%) 1465 (2.7%) 1401 (1.8%) 1596 (11.8%) 1427 

1482 (18.4%) 1326 (6.0%) 1377 (10.1%) 1380 (10.3%) 1265 (1.1%) 1217 (2.7%) 1407 (12.5%) 1251 

1290 (17.2%) 1184 (7.5%) 1208 (9.8%) 1208 (9.7%) 1154 (4.8%) 1102 (0.1%) 1261 (14.5%) 1101 

886 (14.3%) 828 (6.9%) 843 (8.6%) 840 (8.4%) 801 (3.4%) 788 (1.7%) 924 (19.2%) 775 

836 (17.6%) 761 (7.1%) 821 (15.4%) 794 (11.7%) 691 (2.8%) 708 (0.4%) 771 (8.4%) 711 

509 (17.5%) 475 (9.7%) 496 (14.5%) 482 (11.4%) 440 (1.6%) 439 (1.4%) 447 (3.2%) 433 

302 (10.7%) 295 (8.0%) 296 (8.3%) 289 (5.9%) 258 (5.5%) 269 (1.5%) 263 (3.7%) 273 

t lg 1505 (10.9%) 1371 (0.9%) 1403 (3.3%) 1404 (3.4%) 1398 (2.9%) 1346 (0.9%) 1410 (3.8%) 1358 

966 (1.0%) 940 (3.7%) 939 (3.8%) 916 (6.1%) 975 (0.1%) 981 (0.5%) 865 (11.4%) 976 

687 (36.9%) 618 (23.1%) 670 (33.4%) 640 (27.5%) 597 (18.9%) 501 (0.2%) 627 (24.9%) 502 

t 2g 1619 (19.1%) 1453 (6.9%) 1504 (10.6%) 1511 (11.1%) 1470 (8.1%) 1351 (0.7%) 1483 (9.0%) 1360 

912 (0.2%) 907 (0.8%) 898 (1.8%) 868 (5.0%) 890 (2.6%) 931 (1.8%) 919 (0.5%) 914 

903 (4.3%) 701 (19.0%) 828 (4.3%) 734 (15.1%) 834 (3.6%) 847 (2.1%) 784 (9.4%) 865 

647 (14.1%) 637 (12.3%) 634 (11.7%) 613 (8.1%) 637 (12.3%) 541 (4.6%) 591 (4.2%) 567 

Odd Parity 

a u 1111 (2.8%) 1112 (2.5%) 1097 (4.1%) 1061 (7.2%) 1206 (5.5%) 1142 (0.1%) 972 (15.0%) 1143 

g u 1701 (17.7%) 1562 (8.0%) 1607 (11.1%) 1597 (10.4%) 1546 (6.9%) 1413 (2.3%) 1587 (9.8%) 1446 

1527 (16.6%) 1406 (7.3%) 1447 (10.5%) 1437 (9.7%) 1401 (6.9%) 1327 (1.3%) 1436 (9.6%) 1310 

1113 (14.7%) 1031 (6.3%) 1052 (8.4%) 1050 (8.3%) 1007 (3.8%) 961 (0.9%) 1110 (4.1%) 970 

922 (0.2%) 859 (7.0%) 864 (6.5%) 826 (10.6%) 832 (10.0%) 929 (0.5%) 914 (1.1%) 924 

863 (13.5%) 738 (2.9%) 846 (11.4%) 786 (3.4%) 816 (7.4%) 789 (3.8%) 750 (1.3%) 760 

414 (3.5%) 390 (2.4%) 400 (0.1%) 390 (2.4%) 358 (10.5%) 385 (3.8%) 362 (9.5%) 400 

h u 1905 (22.2%) 1762 (13.0%) 1791 (14.9%) 1784 (14.5%) 1646 (5.6%) 1552 (0.4%) 1709 (9.6%) 1559 

1597 (15.3%) 1447 (4.5%) 1487 (7.4%) 1491 (7.7%) 1469 (6.1%) 1385 (0.0%) 1467 (5.9%) 1385 

1453 (30.0%) 1320 (18.2%) 1353 (21.2%) 1354 (21.2%) 1269 (13.6%) 1129 (1.1%) 1344 (20.3%) 1117 

886 (10.5%) 793 (1.0%) 858 (7.1%) 824 (2.8%) 812 (1.4%) 801 (0.0%) 822 (2.6%) 801 

777 (11.7%) 753 (8.2%) 761 (9.3%) 738 (6.0%) 724 (4.0%) 700 (0.6%) 706 (1.4%) 696 

640 (13.6%) 587 (4.2%) 613 (8.9%) 592 (5.2%) 531 (5.7%) 543 (3.6%) 546 (3.0%) 563 

463 (35.0%) 455 (32.7%) 454 (32.4%) 442 (28.7%) 403 (17.5%) 361 (5.0%) 403 (17.5%) 343 

t lu 1637 (14.5%) 1553 (8.6%) 1587 (11.0%) 1549 (8.4%) 1437 (0.6%) 1450 (1.5%) 1628 (13.9%) 1429 

1396 (18.0%) 1245 (5.2%) 1287 (8.8%) 1297 (9.6%) 1212 (2.5%) 1208 (2.1%) 1353 (14.4%) 1183 

656 (13.9%) 614 (6.5%) 623 (8.1%) 625 (8.5%) 637 (10.6%) 589 (2.3%) 719 (24.8%) 576 

627 (18.9%) 575 (9.1%) 621 (17.8%) 599 (13.6%) 544 (3.2%) 505 (4.2%) 577 (9.5%) 527 

t 2u 1828 (15.9%) 1709 (8.4%) 1728 (9.5%) 1713 (8.6%) 1558 (1.2%) 1575 (0.1%) 1687 (7.0%) 1577 
1327 (29.3%) 1214 (1.1%) 1259 (4.5%) 1257 (4.7%) 1241 (3.3%) 1212 (0.9%) 1314 (9.4%) 1201 
1074 (4.7%) 995 (3.1%) 1025 (0.1%) 1014 (1.2%) 999 (2.6%) 1025 (0.1%) 1134 (10.5%) 1026 
835 (22.8%) 765 (12.6%) 830 (22.1%) 799 (17.5%) 690 (1.5%) 677 (0.4%) 776 (14.1%) 680 
393 (10.5%) 377 (5.9%) 385 (8.0%) 372 (4.5%) 350 (1.7%) 367 (3.1%) 348 (2.2%) 356 
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Table VIII: Vibrational frequencies v (cm -1 ) of Ceo calculated by using the B3LYP method with a variety of Pople-style 
basis sets. Numbers in the parenthesis are the relative errors of the calculated frequencies to the experimental frequencies listed 
in Table VII. The approximated CPU times for STO-3G, 3-21G, 6-31G and 6-31G* basis sets are 4 hours, 5 days, 8 days and 
12 days, respectively. The other theoretical results are from Bohnen et al. [87], Dixon et al. [100], Hara et al. [101], and Onida 
et al. [102]. 



Mode 



STO-3G 



3-21G 



6-31G 



6-31G* 



Hara 



Dixon 



Onida 



Bohnen 



Even Parity 
a g 

9g 



Ug 



tig 



1549 (5.4%) 
502 (0.9%) 
1594 (4.5% 
1380 (1.8% 
1127 (4.7% 
788 (2.2% 
592 (4.6% 
508 (4.6% 
1677 (6.3% 
1500 (5.1% 
1332 (6.5% 
1166 (5.9% 
802 (3.5% 
734 (3.3% 
449 (3.7% 
271 (0.6% 
1357 (0.1% 
865 (11.3% 
594 (18.3% 
1431 (5.2% 
827 (9.5% 
809 (6.4% 
581 (2.5% 



Odd Parity 



1501 (2.1%) 
491 (1.4%) 
1524 (0.1%) 
1323 (2.4%) 
1062 (1.3%) 
690 (14.4%) 
598 (3.7%) 
484 (0.4%) 
1609 (2.0%) 
1436 (0.6%) 
1231 (1.6%) 
1112 (1.0%) 
781 (0.8%) 
678 (4.7%) 
429 (1.0%) 
271 (0.8%) 
1278 (5.9%) 
865 (11.4%) 
544 (8.3%) 
1330 (2.2%) 
839 (8.2%) 
650 (24.9%) 
586 (3.3%) 



1524 (3.7%) 
496 (0.4%) 
1546 (1.4%) 
1342 (1.0%) 
1099 (2.1%) 
777 (3.7%) 
594 (4.3%) 
500 (2.9%) 
1627 (3.1%) 
1466 (2.8%) 
1274 (1.8%) 
1129 (2.6%) 
788 (1.7%) 
738 (3.8%) 
448 (3.5%) 
272 (0.5%) 
1303 (4.0%) 
858 (12.1%) 
593 (18.1%) 
1372 (0.8%) 
829 (9.3%) 
766 (11.4%) 
582 (2.6%) 



1504 (2.3%) 
489 (1.8%) 
1538 (0.9%) 
1334 (1.6%) 
1093 (1.6%) 
754 (6.5%) 
577 (7.1%) 
489 (0.7%) 
1618 (2.5%) 
1455 (1.9%) 
1275 (1.9%) 
1125 (2.2%) 
766 (1.2%) 
718 (0.9%) 
436 (0.8%) 
266 (2.7%) 
1301 (4.2%) 
840 (13.9%) 
576 (14.8%) 
1370 (0.7%) 
804 (12.0%) 
743 (14.1%) 
566 (0.1%) 



1531 (4.1%) 
502 (0.8%) 
1538 (0.9%) 
1337 (1.4%) 
1123 4.4(%) 
759 (5.8%) 
579 (6.8%) 
486 (0.0%) 
1609 (2.0%) 
1475 (3.4%) 
1288 (3.0%) 
1129 (2.5%) 
794 (2.5%) 
711 (0.0%) 
430 (0.7%) 
269 (1.5%) 
1305 (3.9%) 
842 (13.7%) 

565 (12.5%) 
1370 (0.7%) 
809 (11.5%) 
765 (11.6%) 

566 (0.2%) 



1525 (3.7%) 
499 (0.2%) 
1548 (1.5% 
1347 (0.7% 
1122 (4.3% 
788 (2.2% 
573 (7.7% 
484 (0.4% 
1618 (2.5% 
1475 (3.4% 
1297 (3.7% 
1128 (2.5% 
788 (1.7% 
727 (2.3% 
431 (0.5% 
261 (4.4% 
1318 (2.9% 
830 (15.0% 
579 (15.3% 
1393 (2.4% 
839 (8.2% 
804 (7.1% 
551 (2.8% 



1447 (1.6%) 
482 (3.2%) 
1479 (3.0%) 
1314 (3.1%) 
1047 (2.7%) 
781 (3.1%) 
594 (4.3%) 
482 (0.8%) 
1573 (0.3%) 
1394 (2.3%) 
1208 (3.4%) 
1098 (0.3%) 
775 (0.0%) 
730 (2.7%) 
435 (0.5%) 
261 (4.4%) 
1284 (5.4%) 
847 (13.2%) 
580 (15.5%) 
1257 (7.6%) 
816 (10.7%) 
789 (8.8%) 
559 (1.4%) 



1475 (0.3%) 
481 (3.4%) 
1501 (1.6%) 
1287 (5.1%) 
1037 (3.6%) 
772 (4.2%) 
570 (8.2%) 
480 (1.2%) 
1580 (0.1%) 
1422 (0.4%) 
1198 (4.2%) 
1079 (2.0%) 
763 (1.5%) 
716 (0.7%) 
422 (2.5%) 
263 (3.7%) 
1241 (8.6%) 
826 (15.4%) 
563 (10.8%) 
1277 (6.1%) 
800 (12.5%) 
788 (8.9%) 
543 (4.4%) 



tii 



t2v 



994 (13.0% 
1513 (4.7% 
1381 (5.4% 
1009 (4.1% 
815 (11.8% 
784 (3.1% 
368 (8.1% 
1668 (7.0% 
1427 (3.0% 
1295 (16.0% 
771 (3.7% 
696 (0.1% 
564 (0.1% 
417 (21.5% 
1505 (5.3% 
1266 (7.0% 
596 (3.4% 
546 (3.6% 
1622 (2.8% 
1235 (2.8% 
1002 (2.4% 
734 (7.9% 
354 (0.5% 



1013 (11.4%) 
1435 (0.8%) 
1303 (0.5%) 
972 (0.2%) 
795 (14.0%) 
666 (12.4%) 
359 (10.3%) 
1596 (2.4%) 
1340 (3.2%) 
1219 (9.1%) 
719 (10.3%) 
674 (3.2%) 
532 (5.5%) 
419 (22.1%) 
1454 (1.7%) 
1175 (0.6%) 
582 (0.9%) 
508 (3.6%) 
1568 (0.6%) 
1163 (3.2%) 
964 (6.0%) 
678 (0.3%) 
345 (3.1%) 



991 (13.3%) 
1470 (1.7%) 
1338 (2.1%) 
985 (1.5%) 
788 (14.7%) 
775 (2.0%) 
366 (8.4%) 
1617 (3.7%) 
1371 (1.0%) 
1246 (11.6%) 
763 (4.7%) 
697 (0.1%) 
555 (1.4%) 
418 (21.8%) 
1479 (3.5%) 
1209 (2.2%) 
587 (1.8%) 
553 (4.9%) 
1579 (0.2%) 
1201 (0.0%) 
984 (4.1%) 
742 (9.2%) 
352 (4.3%) 



982 (14.1%) 
1461 (1.0%) 
1333 (1.7%) 
977 (0.7%) 
787 (14.8%) 
751 (1.2%) 
357 (10.7%) 
1608 (3.2%) 
1369 (11.8%) 
1243 (11.3%) 
736 (8.1%) 
679 (2.5%) 
540 (4.1%) 
408 (19.0%) 
1464 (2.4%) 
1212 (2.5%) 
570 (1.0%) 
537 (1.8%) 
1568 (0.6%) 
1199 (0.2%) 
966 (5.9%) 
722 (6.1%) 
342 (4.0%) 



968 (15.3%) 
1474 (1.9%) 
1345 (2.7%) 
989 (2.0%) 
780 (15.6%) 
762 (0.3%) 
355 (11.3%) 
1598 (2.5%) 
1371 (1.0%) 
1243 (11.3%) 
742 (7.4%) 
677 (2.7%) 
537 (4.6%) 
411 (19.8%) 
1489 (4.2%) 
1222 (3.3%) 
595 (3.3%) 
528 (0.2%) 
1568 (0.6%) 
1231 (2.5%) 
997 (2.8%) 
715 (5.1%) 
343 (3.7%) 



972 (15.0% 
1480 (2.4% 
1359 (3.7% 
984 (1.4%) 
830 (10.2% 
762 (0.3%) 
350 (12.5% 
1611 (3.3% 
1389 (0.3% 
1248 (11.7% 
762 (4.9% 
671 (3.6% 
541 (3.9% 
401 (16.9% 
1486 (4.0% 
1224 (3.5% 
591 (2.6% 
535 (1.5% 
1571 (0.4% 
1234 (2.7% 
996 (2.9% 
726 (6.8% 
342 (3.9% 



934 (18.3%) 
1395 (3.5%) 
1289 (1.6%) 
939 (3.2%) 
796 (13.9%) 
763 (0.4%) 
352 (12.0%) 
1545 (0.9%) 
1314 (5.1%) 
1198 (7.3%) 
769 (4.0%) 
672 (3.4%) 

540 (4.1%) 
404 (17.8%) 
1399 (2.1%) 
1158 (2.1%) 
566 (1.7%) 

541 (2.7%) 
1537 (2.5%) 
1108 (7.7%) 
936 (8.8%) 
774 (13.8%) 
340 (4.5%) 



973 (14.9%) 
1420 (1.8%) 
1259 (3.9%) 
937 (3.4%) 
790 (14.5%) 
756 (0.5%) 
348 (13.0%) 
1566 (0.4%) 
1291 (6.8%) 
1175 (5.2%) 
750 (6.4%) 
661 (5.0%) 
527 (6.4%) 
388 (13.1%) 
1457 (2.0%) 
1143 (3.4%) 
569 (1.2%) 
514 (2.5%) 
1546 (2.0%) 
1131 (5.8%) 
945 (7.9%) 
725 (6.6%) 
343 (3.8 
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Table IX: Vibrational frequencies v (cm -1 ) of Ceo obtained by Choi et al. [141] using B3LYP/3-21G but involving scaling 
of the internal force constant by using Pulay's method. Numbers in the parenthesis are the relative errors of the calculated 
frequencies to the experimental frequencies listed in Table VII. 







Even Parity 










Odd Parity 






Mode 


V 


Mode v 


Mode 


V 


Mode 


V 


Mode v 


Mode 


V 


a g 


1470 (0.0%) 




hg 


1290 (5.0%) 


a u 


1078 (5.7%) 




hu 


1433 (0.3%) 




495 ( 0.6%) 


h g 1576 (0.1%) 




904 (7.4%) 










1180 (0.3%) 






1427 (0.0%) 




565 (12.5%) 






h u 1567 (0.5%) 




577 (0.2%) 


9g 


1497 (1.8%) 


1251 (0.0%) 






9u 


1429 (1.2%) 


1343 (3.0%) 




526 (0.2%) 




1348 (0.6%) 


1101 (0.0%) 








1315 (0.4%) 


1214 (8.7%) 


hu 


1524 (3.4%) 




1040 (3.3%) 


775 (0.0%) 




1340 (1.5%) 




970 (0.0%) 


737 (8.0%) 




1142 (5.0%) 




758 (6.0%) 


711 (0.0%) 




831 (9.1%) 




797 (13.7%) 


694 (0.3%) 




955 (6.9%) 




592 ( 4.7%) 


431 (0.5%) 




668 (22.8%) 




707 (7.0%) 


535 (5.0%) 




716 (5.3%) 




485 ( 0.2%) 


267 (2.2%) 




614 (8.3%) 




354 (11.5%) 


403 (17.5%) 




340 (4.5%) 



Table X: Fifty eight IR-active frequencies {y, in cm x ) of C4SN12 calculated by using RHF and B3LYP methods with a 
variety of Pople-style basis sets. 







Doubly-degenerate Modes 










Non-degenerate Modes 






B3LYP 






RHF 


B3LYP 






RHF 


STO-3G 








1 U-oO 


o-ZLKj 




6-31G* 


STO-3G 


o-ZLKj 








i-ZLKj 




6-31G* 


314 


308 


318 


309 


O K O 


338 


350 


344 


286 


287 


297 


007 
lot 


OAK 

305 


303 


317 


305 


363 


359 


369 


357 


400 


386 


398 


381 


328 


318 


329 


325 


370 


350 


362 


361 


394 


404 


406 


395 


430 


431 


436 


419 


362 


359 


371 


360 


404 


394 


406 


392 


406 


412 


418 


406 


452 


452 


454 


443 


391 


402 


406 


395 


428 


434 


437 


423 


437 


422 


441 


425 


513 


472 


495 


479 


450 


413 


447 


431 


527 


466 


498 


482 


465 


440 


464 


446 


551 


497 


523 


506 


470 


454 


467 


460 


556 


505 


524 


522 


496 


469 


491 


481 


578 


532 


554 


549 


581 


567 


580 


579 


641 


603 


617 


614 


566 


558 


569 


570 


619 


591 


606 


604 


592 


600 


639 


611 


695 


662 


720 


662 


623 


585 


633 


606 


723 


662 


707 


671 


645 


609 


653 


622 


735 


709 


739 


705 


632 


635 


662 


639 


738 


730 


741 


719 


656 


654 


674 


645 


767 


734 


747 


722 


656 


665 


679 


658 


763 


747 


755 


735 


692 


698 


707 


686 


802 


762 


779 


755 


675 


686 


697 


675 


771 


767 


786 


757 


735 


709 


722 


700 


839 


800 


819 


770 


695 


701 


709 


695 


801 


783 


801 


779 


784 


798 


798 


778 


871 


859 


867 


836 


722 


721 


729 


712 


843 


810 


829 


807 


798 


809 


803 


786 


899 


878 


879 


853 


788 


804 


800 


782 


879 


871 


874 


847 


958 


963 


968 


961 


1058 


1035 


1044 


1030 


973 


958 


965 


954 


1072 


1028 


1040 


1020 


972 


977 


982 


967 


1076 


1055 


1062 


1045 


1055 


1041 


1050 


1040 


1162 


1118 


1135 


1116 


990 


1010 


993 


976 


1095 


1102 


1092 


1060 


1203 


1136 


1163 


1159 


1336 


1217 


1259 


1244 


1057 


1042 


1056 


1047 


1165 


1125 


1145 


1124 


1231 


1172 


1199 


1197 


1372 


1265 


1301 


1292 


1220 


1142 


1171 


1177 


1351 


1232 


1272 


1272 


1275 


1210 


1237 


1237 


1415 


1304 


1341 


1334 


1276 


1215 


1238 


1229 


1393 


1311 


1345 


1330 


1348 


1280 


1306 


1302 


1501 


1369 


1409 


1397 


1331 


1252 


1277 


1270 


1446 


1326 


1365 


1351 


1391 


1311 


1339 


1340 


1552 


1413 


1454 


1450 


1362 


1293 


1318 


1317 


1464 


1381 


1421 


1377 


1440 


1354 


1387 


1389 


1618 


1467 


1512 


1511 


1421 


1345 


1380 


1384 


1521 


1405 


1463 


1412 


1478 


1410 


1448 


1437 


1647 


1511 


1562 


1540 


1437 


1374 


1413 


1394 


1609 


1451 


1502 


1494 


1525 


1443 


1476 


1473 


1731 


1566 


1613 


1605 


1466 


1381 


1422 


1418 


1621 


1490 


1532 


1533 


1555 


1475 


1502 


1499 


1769 


1612 


1652 


1642 


1506 


1418 


1458 


1461 


1700 


1523 


1583 


1574 


1602 


1521 


1553 


1553 


1839 


1673 


1713 


1713 


1560 


1497 


1529 


1518 


1773 


1633 


1674 


1662 


1628 


1544 


1578 


1576 


1855 


1693 


1735 


1748 


1617 


1516 


1549 


1555 


1872 


1692 


1731 


1748 


1680 


1592 


1625 


1624 


1941 


1760 


1801 


1808 


1627 


1542 


1574 


1575 


1894 


1701 


1742 


1751 
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Table XI: Fifty eight Raman-active frequencies (y, in cm 1 ) for C48N12 calculated by using RHF and B3LYP methods with 
a variety of Pople-style basis sets. 



Doubly-degenerate Modes 

B3LYP RHF ~ 

ST0-3G 3-21G 6-31G 6-31G* STO-3G 3-21G 6-31G 6-31G* 



245 


248 


252 


245 


263 


262 


268 


261 


259 


261 


265 


260 


286 


281 


286 


281 


376 


368 


387 


371 


429 


409 


428 


413 


410 


388 


409 


396 


470 


433 


455 


440 


447 


429 


446 


437 


514 


479 


497 


488 


493 


450 


482 


472 


583 


517 


549 


540 


549 


568 


566 


551 


604 


610 


611 


592 


580 


596 


597 


581 


652 


647 


650 


629 


617 


613 


652 


627 


725 


699 


736 


706 


649 


629 


662 


641 


773 


726 


759 


741 


665 


671 


698 


671 


780 


762 


785 


757 


716 


695 


714 


699 


831 


776 


807 


788 


767 


760 


768 


766 


848 


820 


833 


821 


782 


776 


782 


780 


871 


844 


853 


840 


841 


865 


860 


843 


935 


943 


944 


918 


858 


887 


877 


854 


961 


974 


966 


933 


1124 


1084 


1105 


1093 


1238 


1160 


1192 


1164 


1162 


1123 


1141 


1138 


1281 


1197 


1225 


1215 


1182 


1147 


1163 


1156 


1313 


1238 


1262 


1250 


1274 


1186 


1220 


1223 


1423 


1276 


1325 


1318 


1318 


1243 


1274 


1272 


1476 


1341 


1383 


1377 


1386 


1311 


1337 


1336 


1546 


1421 


1457 


1449 


1452 


1369 


1407 


1404 


1608 


1468 


1517 


1503 


1466 


1388 


1420 


1420 


1647 


1490 


1537 


1540 


1504 


1415 


1452 


1451 


1717 


1541 


1591 


1587 


1571 


1489 


1517 


1515 


1781 


1631 


1671 


1664 


1595 


1516 


1551 


1551 


1831 


1670 


1710 


1720 


1633 


1544 


1578 


1578 


1883 


1708 


1749 


1759 


1662 


1578 


1610 


1610 


1916 


1735 


1778 


1788 



Non-degenerate Modes 

B3LYP RHF ~ 

STO-3G 3-21G 6-31G 6-31G* STO-3G 3-21G 6-31G 6-31G* 



264 


264 


268 


264 


294 


288 


291 


288 


382 


368 


389 


376 


445 


406 


429 


388 


406 


398 


415 


398 


463 


415 


448 


425 


444 


415 


442 


424 


510 


455 


483 


454 


471 


458 


477 


467 


523 


491 


511 


489 


498 


491 


500 


495 


576 


521 


543 


524 


544 


505 


544 


510 


608 


566 


587 


575 


578 


588 


592 


576 


646 


637 


640 


615 


585 


598 


602 


588 


650 


651 


652 


637 


621 


603 


616 


597 


723 


655 


670 


652 


625 


614 


652 


627 


744 


705 


749 


716 


658 


657 


684 


656 


777 


760 


784 


738 


730 


690 


709 


689 


832 


779 


802 


773 


768 


765 


773 


766 


867 


825 


836 


819 


836 


851 


844 


830 


924 


918 


919 


897 


856 


881 


869 


844 


953 


965 


956 


916 


1115 


1080 


1100 


1092 


1216 


1165 


1192 


1175 


1194 


1160 


1173 


1162 


1328 


1244 


1264 


1251 


1221 


1181 


1196 


1189 


1344 


1260 


1281 


1273 


1316 


1241 


1267 


1265 


1400 


1335 


1374 


1348 


1337 


1253 


1280 


1279 


1469 


1347 


1389 


1377 


1367 


1293 


1322 


1320 


1498 


1375 


1424 


1380 


1416 


1347 


1383 


1379 


1534 


1415 


1468 


1429 


1439 


1373 


1415 


1397 


1583 


1437 


1485 


1476 


1481 


1419 


1453 


1441 


1655 


1528 


1578 


1552 


1531 


1440 


1477 


1477 


1721 


1560 


1612 


1602 


1558 


1487 


1511 


1505 


1766 


1612 


1649 


1626 


1581 


1509 


1545 


1530 


1817 


1644 


1687 


1680 


1683 


1592 


1624 


1623 


1946 


1767 


1804 


1809 



Table XII: RHF and B3LYP calculations of IR intensities (Iir, in 10 3 m/mole ) of C60 with the corresponding vibrational 
modes and frequencies ( v, in cm -1 ). 



Method 


Mode 


STO-3G 


3-21G 




6-31G 




6-31G 


•f 


Iir 


V 


Iir 


V 


Iir 


V 


Iir 


V 


B3LYP 


tlu 


10.6 


1505 


14.0 


1454 


17.2 


1479 


15.6 


1464 






21.0 


1266 


9.2 


1175 


10.1 


1209 


8.9 


1212 






0.5 


596 


5.9 


582 


8.2 


587 


10.7 


570 






22.0 


546 


28.8 


508 


27.7 


553 


27.1 


537 


RHF 


tlu 


11.5 


1637 


16.5 


1553 


24.5 


1587 


17.1 


1549 






25.1 


1396 


11.4 


1245 


11.8 


1287 


10.6 


1297 






0.6 


656 


6.0 


614 


14.7 


623 


10.4 


625 






35.9 


627 


42.1 


575 


34.3 


621 


48.0 


599 
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FIG. 4: Ab initio calculation of IR intensity (Iir, in 10 3 m/mole) at its corresponding frequency in C48N12: 
(a) B3LYP/STO-3G; (b) B3LYP/3-21G; (c) RHF/3-21G. The solid lines are the calculated results for C 60 . 



shows collective vibration along the z-x direction, while 
the pentagon structure for the high-frequency case con- 
tracts, accompanying a large stretching of site 3 along 
the z direction. 



AZ 



4 




FIG. 5: The vibrational displacements of 
sites 1 to 5 for the strongest IR spectral sig- 
nals of C48N12 in the low-frequency (filled 
circles) and high-frequency (open circles) re- 
gions for the B3LYP/3-21G case. 



VI. NUCLEAR MAGNETIC SHIELDING 
TENSORS 

A. Theoretical Methods 

There are a number of theoretical methods for calcu- 
lating the second-order magnetic response properties of 
molecules. It has been shown that gauge-invariant proce- 
dures (GIPs) [71,143] are required to predict accurately 
these properties within a finite basis approximation. In 



this paper, we focus on using two GIPs, i.e., the gauge- 
including atomic orbital (GIAO) procedure and the con- 
tinuous set of gauge transformations (CSGT) procedure, 
to predict NMR shielding tensors at the Hartree-Fock 
and DFT levels of theory GIAO and CSGT achieve 
gauge-invariance in different ways. The GIAO method 
uses basis functions having explicit field dependence [70] , 
whereas the CSGT method achieves the gauge-invariance 
by performing a continuous set of gauge transformations. 
Ditchfield [144] first adopted the GIAO method to per- 
form quantum chemical NMR shift calculations and the 
CSGT method was developed by Keith and Bader [71]. 
In the following, we briefly introduce the two methods. 
More details can be found, for example, in the works of 
Cheeseman et al. [145] and Wolinski et al. [146]. 

It is known that the nuclear magnetic shielding tensor 
can be written as the mixed second derivative of the en- 
ergy E with respect to the external magnetic field B and 
the magnetic moment I of nucleus X: 

where Bj, and 1^ are the components of the external mag- 
netic field and induced magnetic moment, respectively. 
The nuclear magnetic shielding isotropy o is defined as 
[65] 

V = {?xx + Oyy + <?zz)/3, (12) 

and the shielding anisotropy Act, an indication of the 
quality of the magnetic shielding tensor, is defined as 
[65] 

A<7 = <7 3 - ((7! + (T 2 )/2, (13) 

where u\ < a 2 < 03 are the eigenvalues of the sym- 
metrized shielding tensor. The nuclear magnetic shield- 
ing difference, or say, chemical shift 5, is reported in ppm 
(i.e., parts per million) and given by [65] 

g _ ^(reference) _ a (sample)\ x ^q6 ^ 
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Where preference) and ^sample) denote tne shielding 

isotropics a for the reference and sample, respectively. 

In both DFT and HF theory, the expression for the 
nuclear magnetic shielding tensor of nucleus X is given 
by [145] 

dO 

a lk =< U ik O > + < k -— > (15) 

with 

tt /* , r * ( r - ^x)kk - n(r - R x )k , - . , . 
n, fe = ($J 2c 2|r-R x |3 (16) 



i|r-Rx) x y]fc 

where c is the velocity of light, O is the density matrix, 
<I> M and are spin orbitals, Hx is the position vec- 
tor of nuclear X, and r is the real space vector. In the 
above equation, dO/dB k is the derivative of the density 
matrix O with respect to the fcth component of the mag- 
netic field B and is obtained via solution of the coupled- 
perturbed equations [145] for the appropriate perturba- 
tion. As gauge-invariance is achieved in different ways, 
the GIAO and CSGT methods differ at this point in the 
formation of the coupled-perturbed equations [145]. 



Table XIII: RHF and B3LYP calculations of the absolute isotropy (a, in ppm) and anisotropy (Act, in ppm) of the carbon 
and nitrogen shielding tensors with a variety of Pople-style basis sets for C48N12 aza-fullcrenc, Ceo and tetramethylsilane (TMS) 
by using the GIAO method. Numbers in the parenthesis for Ceo are the relative errors of the calculated 13 C NMR shift 8 to 
the NMR chemical shift 5 exp - = 142.7 ppm measured by Taylor et aJ. [149]. 











STO-3G 


3-21G 


6-31G 


6-31G* 


Method 


Molecule 


Site Numbers {rii} 


Nuclei 


a Aa 


a Aa 


a Aa 


a Aa 



RHF C48N12 {1, 13, 16, 31, 38, 51} 9L5 155.5 65T9 155.6 441) 173.0 48J3 168.3 

{2, 12, 29, 32, 37, 52} 13 C 109.8 156.7 89.8 149.2 67.9 169.2 67.7 168.2 

{3, 11, 28, 33, 36, 53} 13 C 116.8 132.2 99.2 123.5 76.6 141.2 88.5 126.9 

{4, 15, 27, 34, 40, 54} 13 C 103.5 141.9 75.8 139.0 57.9 154.9 46.9 156.7 

{ 5, 14, 30, 35, 39, 55} 15 N 193.7 134.8 142.7 143.8 114.4 171.5 140.7 128.5 

{6, 18, 24, 42, 48, 58} 13 C 120.2 104.8 101.6 89.8 78.2 110.5 89.6 82.5 

{7, 19, 23, 43, 47, 57} 13 C 106.2 135.9 80.6 131.0 59.4 152.2 54.2 148.7 

{8, 20, 22, 44, 46, 56} 13 C 111.5 140.9 94.6 126.1 70.7 150.8 82.9 131.1 

{9, 21, 26, 45, 50, 60} 15 N 182.7 151.6 132.6 167.7 95.3 202.6 121.7 170.4 

{10, 17, 25, 41, 49, 59} 13 C 104.9 125.2 80.0 117.0 59.5 133.5 55.1 126.9 

C 60 {1, 2, 3, 4, 5, 6, 60} 13 C 101.8 156.6 74.0 162.7 54.1 180.1 54.7 178.8 

TMS carbon site 13 C 238.5 5.6 208.3 16.9 200.8 21.4 195.1 17.5 

NH 3 nitrogen site 15 N 306.9 9.2 271.1 17.5 264.4 18.2 260.8 17.4 

Ceo Calculated 13 C NMR shift S 136.7 (4.2%) 134.3 (5.9%) 146.7 (2.8%) 140.4 (1.6%) 



B3LYP 



C48N12 


{1, 13, 16, 31, 38, 51} 


13 C 


103.9 


109.5 


70.9 


117.5 


53.3 


128.0 


51.9 


130.4 




{2, 12, 29, 32, 37, 52} 


13 C 


114.6 


121.9 


84.6 


123.6 


66.4 


138.1 


62.8 


138.9 




{3, 11, 28, 33, 36, 53} 


13 C 


120.6 


101.1 


91.8 


98.6 


74.5 


110.9 


72.9 


107.3 




{4, 15, 27, 34, 40, 54} 


13 C 


113.3 


103.6 


82.8 


104.6 


66.3 


117.1 


61.8 


113.9 




{ 5, 14, 30, 35, 39, 55} 


15 N 


173.2 


126.2 


116.4 


143.2 


92.8 


155.4 


103.8 


139.3 




{6, 18, 24, 42, 48, 58} 


13 C 


123.5 


70.4 


92.5 


66.5 


76.1 


76.4 


74.3 


63.0 




{7, 19, 23, 43, 47, 57} 


13 C 


114.1 


96.4 


80.6 


100.2 


63.0 


76.4 


57.4 


114.0 




{8, 20, 22, 44, 46, 56} 


13 c 


119.5 


101.7 


90.5 


104.7 


73.8 


116.2 


72.3 


116.2 




{9, 21, 26, 45, 50, 60} 


15 N 


171.7 


121.9 


110.4 


142.2 


84.0 


157.5 


91.2 


146.7 




{10, 17, 25, 41, 49, 59} 


13 C 


113.9 


121.9 


82.0 


84.6 


65.6 


94.4 


61.7 


90.0 


Ceo 


{1, 2, 3, 4, 5, 6, 60} 


13 C 


106.1 


125.8 


70.2 


139.3 


53.6 


152.2 


50.5 


154.2 


TMS 


carbon site 


13 c 


225.3 


9.6 


197.6 


20.0 


188.3 


25.2 


183.8 


22.3 


NH 3 


nitrogen site 


15 N 


291.3 


7.3 


262.6 


18.9 


256.8 


21.0 


254.9 


19.3 


Ceo 


Calculated 13 C NMR shift 5 


119.2 (16.5%) 


127.4 (10.7%) 


134.7 (5.6%) 


133.3 


(6.6%) 



Experiment 



Ceo 
TMS 
NH 3 



13 C NMR chemical shift ,5 from Taylor et al. [149]: 
absolute shielding isotropy a for 13 C from Jameson et al. [148]: 
absolute shielding isotropy (a) /anisotropy (Aa) for N from Ref. [150,151]: 



142.7 
188.1 
264.5/20 
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1. GIAO Method 



2. CSGT Method 



The GIAO method uses the following explicit field- 
dependent basis functions for calculating the magnetic 
shielding tensor: 



- „-i(Bxr M ).r/(2c) 



*m(0)). 



(18) 



where is the position vector of basis function (f^ and 
$ M (0) denotes the usual field-independent basis func- 
tions. In this method, three sets of the coupled-perturbed 
equations [145] are solved, one for each direction of the 
magnetic field. 



In the CSGT method, the gauge-invariance is achieved 
from accurate calculations of the induced first-order elec- 
tronic current density J^'(r) by performing a gauge 
transformation for each point in space. The nuclear mag- 
netic shielding tensor is expressed in terms of J^^(r), i.e., 



Oik 



1 

Be 



Rxxj["(r) 



dRx, 



(19) 



Table XIV: RHF and B3LYP calculations of the absolute isotropy (ct, in ppm) and anisotropy (Act, in ippm) of the carbon 
and nitrogen shielding tensors with a variety of Pople-style basis sets for C48N12 aza-fullerene, Ceo and tetramethylsilane (TMS) 
by using CSGT methods. Numbers in the parenthesis for C60 are the relative errors of the calculated 13 C NMR shift <5 to the 
NMR chemical shift 5 exp - = 142.7 ppm measured by Taylor et aJ. [149]. 



Method 


Molecule 


Site Numbers {rii} 


Nuclei 


STO-3G 


3-21G 


6-31G 


6-31G* 


CT 


Act 


CT 


Act 


CT 


Act 


CT 


Act 


RHF 


C4SN12 


{1, 13, 16, 31, 38, 51} 


i3 C 


45.3 


125.0 


60.7 


151.5 


30.7 


168.9 


41.3 


170.1 






{2, 12, 29, 32, 37, 52} 


13 C 


55.8 


127.0 


84.1 


143.8 


54.0 


162.1 


60.4 


168.3 






{3, 11, 28, 33, 36, 53} 


13 C 


60.0 


114.0 


91.2 


122.4 


60.9 


139.1 


78.9 


129.8 






{4, 15, 27, 34, 40, 54} 


13 c 


52.9 


119.5 


68.8 


138.7 


41.1 


155.5 


39.2 


160.4 






{ 5, 14, 30, 35, 39, 55} 


13 N 


78.0 


134.0 


120.4 


154.3 


86.3 


176.5 


128.2 


136.4 






{6, 18, 24, 42, 48, 58} 


13 C 


63.8 


97.9 


89.8 


96.5 


59.3 


115.5 


79.8 


87.3 






{7, 19, 23, 43, 47, 57} 


13 C 


55.5 


113.5 


72.2 


130.9 


43.2 


149.5 


46.5 


150.5 






{8, 20, 22, 44, 46, 56} 


13 c 


58.6 


118.0 


84.9 


125.5 


52.4 


150.1 


73.7 


135.2 






{9, 21, 26, 45, 50, 60} 


13 N 


74.0 


144.8 


116.0 


170.8 


72.2 


203.6 


112.2 


173.1 






{10, 17, 25, 41, 49, 59} 


13 C 


55.3 


107.8 


73.1 


118.7 


43.3 


136.9 


47.3 


129.5 




Ceo 


{1,2, 3, 4, 5, 6, ...,60} 


13 C 


50.2 


128.8 


71.2 


156.1 


43.3 


172.7 


48.9 


178.6 




TMS 


carbon site 


13 c 


131.9 


9.5 


191.0 


6.3 


185.5 


13.1 


190.6 


13.7 




NH 3 


nitrogen site 


13 N 


148.1 


18.6 


200.2 


5.1 


209.2 


6.3 


228.4 


6.6 




Ceo 


Calculated 13 C NMR shift 6 


81.7 (42.7%) 


119.8 


(16.0%) 


142.2 


(0.4%) 


141.7 


(0.7%) 


B3LYP 


C48N12 


{1, 13, 16, 31, 38, 51} 


13 C 


58.6 


86.4 


66.0 


114.1 


39.1 


126.6 


45.7 


131.8 






{2, 12, 29, 32, 37, 52} 


13 C 


66.1 


94.8 


81.7 


119.3 


53.9 


133.5 


57.2 


139.6 






{3, 11, 28, 33, 36, 53} 


13 C 


69.3 


83.5 


87.2 


99.0 


60.3 


111.1 


66.3 


109.4 






{4, 15, 27, 34, 40, 54} 


13 c 


64.6 


84.8 


75.8 


104.0 


49.1 


118.5 


54.7 


117.6 






{ 5, 14, 30, 35, 39, 55} 


13 N 


78.5 


116.6 


104.9 


142.0 


74.6 


152.7 


96.8 


142.4 






{6, 18, 24, 42, 48, 58} 


13 C 


71.9 


63.6 


84.1 


72.0 


58.7 


81.9 


66.9 


66.5 






{7, 19, 23, 43, 47, 57} 


13 C 


66.5 


77.7 


73.8 


99.7 


47.5 


112.2 


51.0 


115.3 






{8, 20, 22, 44, 46, 56} 


13 C 


70.4 


81.9 


83.5 


101.7 


56.7 


115.8 


65.7 


118.7 






{9, 21, 26, 45, 50, 60} 


13 N 


81.6 


112.8 


103.4 


140.2 


69.3 


156.6 


86.0 


148.5 






{10, 17, 25, 41, 49, 59} 


13 C 


66.9 


71.6 


76.1 


88.1 


49.5 


100.1 


54.8 


93.8 




Ceo 


{1, 2, 3, 4, 5, 6, 60} 


13 c 


58.2 


100.3 


69.2 


134.7 


43.1 


148.4 


45.9 


149.2 




TMS 


carbon site 


13 c 


124.9 


10.0 


181.7 


9.0 


175.0 


16.5 


181.6 


18.0 




NH 3 


nitrogen site 


13 N 


137.7 


27.1 


190.0 


9.6 


200.4 


6.0 


221.5 


6.5 




Ceo 


Calculated 13 C NMR shift 6 


66.7 (53.3%) 


112.5 


(21.2%) 


131.9 


(7.6%) 


135.7 


(4.9%) 



where B is the magnitude of magnetic field. Since the ba- solved (three for the components of the angular momen- 
sis functions do not depend on the magnetic field, only six turn perturbation using any single gauge origin, and three 
sets of the coupled-perturbed equations [145] are left and for the linear momentum perturbation resulting from any 
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single shift in gauge origin). Furthermore, the shielding 
tensors are obtained via Becke's multi-center numerical 
integration scheme [147] 



B. Results 

The above mentioned GIAO and CSGT methods have 
been implemented into the Gaussian 98 program [73,72] 
and all calculations in this section are performed using 
this program. 

Here we perform RHF and B3LYP hybrid DFT calcu- 
lations of the nuclear magnetic shielding tensor of carbon 
and nitrogen atoms in C48N12, Cgo, and tetramcthylsi- 
lane (TMS [148]) by using GIAO and CSGT methods 
with STO-3G, 3-21G, 6-31G and 6-31G* basis sets. Cal- 
culation performed with a specific basis set and ab ini- 
tio method use geometry optimized with the same basis 
set and ab initio method (see section II). Detailed re- 
sults are summarized in Table XIII and XIV. Since no 
present functional includes a magnetic field dependence, 
the DFT methods do not provide systematically better 
shielding results than RHF [72]. Nevertheless, we also list 
the DFT results in Table XIII and XIV for comparison. 
The isotropic NMR chemical shifts 5 relative to that of 
TMS are also shown in Table XIII and XIV. It is seen that 
the NMR shielding tensors for C48N12 are separated into 
eight groups for carbon atom and two groups for nitrogen 
atoms, i.e., eight 13 C and two 15 N NMR spectral signals 
are predicted. In contrast, Ceo has only one 13 C NMR 
signal which is in agreement with experiment [149]. The 
experimental values of the absolute shielding constant ct 
for 15 N in NH 3 and 13 C in TMS, as shown in Table XIII 
and XIV, are 264.5 ppm [150,151] and 188.1 ppm [148], 
and our best calculated results for those reference mate- 
rials are in good agreement with the experiments. 

Table XIII and XIV also demonstrate the convergence 
of the GIAO and CSGT methods with respect to basis set 
for absolute shielding constants calculated with RHF and 
B3LYP hybrid DFT methods. For each first principles 
method, the shielding constants calculated with GIAO 
and CSGT methods are found to converge to almost the 
same values for the large basis set 6-3 1G*. NMR chemi- 
cal shifts, as opposed to the absolute shielding constants 
ct, are measured with high accuracy in applications of 
NMR spectroscopy. On the other hand, calculated chem- 
ical shifts are in better agreement with experiment as 
relative differences arc better represented [145]. Hence, 
given the absolute shielding constants for 13 C and 15 N 
in reference materials, i.e., TMS and NH 3 shown in Ta- 
ble XIII and XIV, we are able to calculate the chem- 
ical shifts 5 and compare them with the experiments. 
For Ceo, we find that the NMR chemical shifts 5 calcu- 
lated with the CSGT method at the RHF/6-31G and 
RHF/6-31G* levels are in agreement with experiment 

(gexp = 142 y ppm measured by Taylor et al [ 14 g])_ Xhis 

suggests that CSGT method would predict NMR spec- 



tral signals much better than GIAO. However, CSGT 
costs more CPU times than GIAO on the same machine, 
for example, about 5 hours more for RHF/6-31G* and 2 
hours for B3LYP/6-31G* in the CSGT method. 

To predict accurate NMR chemical shifts for large 
molecules, it is necessary to assess the accuracy of the 
available ab initio methods by employing lower levels of 
theory and by using basis sets as small as possible [145]. 
The results in Table XIII and XIV indicate that the RHF 
method yields better chemical shifts than the B3LYP hy- 
brid DFT method, especially, for GIAO method and the 
minimum basis set STO-3G. As shown in Table XIII, the 
13 C chemical shift in Cgo differs from experiment by 23.5 
ppm at the B3LYP level, while it differs by only 6 ppm 
at the RHF level of theory. Also it is noted that RHF 
makes positive and negative errors, while B3LYP makes 
only positive errors. 

In addition, nuclear magnetic shielding anisotropics 
Act are reported in Table XIII and XIV. Although these 
properties usually cannot be determined experimentally 
in the gas phase (except in cases where the high sym- 
metry of the molecule together with the calculated dia- 
magnetic contribution to the shielding tensor allows the 
determination of the anisotropy from the spin rotation 
[152]), these values are of interest. Anisotropics can 
be determined in solid state NMR experiments and cal- 
culations are often important for a correct assignment. 
From Table XIII and XIV, we find that the shielding 
anisotropics Act for 15 N in NH 3 calculated with GIAO 
method with both ab initio theories are in good agree- 
ment with experiment, while those calculated with CSGT 
agree poorly with experiment. 



VII. SUMMARY AND OUTLOOK 

In summary, we have performed ab initio calculations 
of the structures, electronic properties, vibrational fre- 
quencies, IR intensities, NMR shielding tensors, linear 
polarizabilities, first- and second-order hyperpolarizabil- 
ities of the C48N12 aza-fullcrenc. Calculated results of 
C48N12 are compared with those of Ceo at the same level 
of theory. It is found that this aza-fullerene has some 
remarkable features, which are different from and can 
compete with C@o- The detailed studies of Cqq show the 
importance of electron correlations and the choice of ba- 
sis sets in the ab initio calculations. Our best results 
for C 60 obtained with the B3LYP hybrid DFT method 
are in excellent agreement with experiment and demon- 
strate the needed efficiency and accuracy of this method 
for obtaining quantitative information on the structural, 
electronic and vibrational properties of these materials. 

Laser sources are widely used in the laboratory and 
industry. However, they are a potential hazard for eyes, 
thermal cameras and other optical sensors [10,11,153]. 
Development of optical limiters which can suppress un- 
desired radiation and effectively decrease transmittancc 
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at high intensity or fluence is necessary. An ideal optical 
limiter [10,153] should have reasonable linear transmit- 
tance at low input fluence (at least of 70%), protect- 
ing optical sensors or eyes against laser pulses of any 
wavelength and pulse duration, and its output energy 
must remain below the optical damage threshold of sen- 
sors or eyes at high fluences. Ceo and fullercnc deriva- 
tives are good candidates for optical limiting applications 
[10,154,155]. In this paper, we found that the average 
second-order hyperpolarizability of C48N12 is about 55% 
larger than that of C 6 o • Hence, it is expected that C4 8 N 12 
is also a good candidate for optical limiting applications. 

Non-linear optical (NLO) materials have vast techno- 
logical applications in telecommunications, optical data 
storage and optical information processing. The search 
for materials with such properties is the subject of intense 
research [9,10,156]. Organic molecules, which have high 
NLO response, could be designed by linking appropriate 
electron-donor and acceptor groups with a spacer made 
of one or several units of conjugated molecules such as a 
polyene or an aromatic chain [156]. Based on such donor- 
acceptor model, we can link donor C4 8 N 12 and acceptor 
Ceo with a polyene or an aromatic chain and design one 
kind of organic molecules. In such organic molecules, 
charge would migrate from C48N12 to Ceo upon electronic 
excitation, giving rise to a large dipole moment along 
the direction connecting the donor /acceptor pair. Thus, 
large NLO response is expected in these donor/acceptor- 
based molecules. 

Carbon nanotubes are promising materials for build- 
ing electronic devices, in particular, field effect tran- 
sistors (FETs) [11,157]. However, single-walled carbon 
nanotube (SWNT) FETs built from as-grown tubes are 
unipolar p-type, i.e., there is no electron current flow even 
at large positive gate biases. This behavior suggests the 
presence of a Schottky barrier at the metal-nanotube con- 
tact. Obviously, the capability to produce n-type tran- 
sistors is important technologically, as it allows the fab- 
rication of nanotube-based complementary logic devices 
and circuits [158,159]. Very recently, Bockrath et al. [160] 
have reported a controlled chemical doping of individual 
semiconducting SWNT ropes by reversibly intercalating 
and deintercalating potassium. Potassium doping could 
change the carriers in the ropes from holes to electrons 
[160]. Their experiments open the way toward other ex- 
periments that require controlled doping such as making 
nanoscale p-n junctions. As shown in this paper, C48N12 
is a good electron donor. We found that incorporating 
C 48 Ni2 into a (10,10) semiconducting SWNT [7] would 
result in -0.26 e charge [90] on the nanotube forming a 
n-type SWNT-based transistor. 

To obtain molecular rectification, the LUMO of the ac- 
ceptor should lie at or slightly above the Fermi level of the 
electrode and above the HOMO of the donor [161]. Hence 
it is important to search for desired donor/acceptor pairs 
which satisfy those conditions. The acceptor/donor pair, 
C48B12/C48N12, actually is demonstrated to be an ideal 
candidate for molecular electronics [90]. For example, 



we have shown that a rectifier formed by C48B12/C48N12 
pair shows a perfect rectification behavior [90]. 

As shown before, C48B12 and C48N12 fee solids are 
semiconducting materials. Because both C48B12 and 
C48N12 molecules have opposite electronic polarizations, 
it is possible to use them to build semiconducting mate- 
rials with opposite electronic polarizations. 

Therefore, C48N12 could have potential applications as 
semiconductor components, nonlinear optical materials, 
and possible building materials for molecular electron- 
ics and photonic devices. Efficiently synthesizing C48N12 
would be of great experimental interest within reach of 
today's technology. 
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